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THE  PRECESSIONAL  CONSTANT. 


To  the  first  of  these  is  frequently  applied  the  term  "Constant  of  Precession." 
But  as  it  is  a  variable  quantity,  the  application  of  the  term  "  constant "  to  it  might  be 
misleading.  Indeed,  no  one  of  the  four  quantities  just  described  can  be  regarded  as 
a  fundamental  constant  in  astronomy,  not  only  because  they  all  vary,  but  because  at 
least  two  of  them  are  a  result  of  both  motions — that  of  the  ecliptic  and  that  of  the 
equator. 

In  the  case  of  the  lunisolar  precession,  defined  as  the  movement  of  the  equinox 
upon  a  fixed  ecliptic,  there  is  some  ambiguity.  As  the  term  is  used  by  Bessel  and 
others,  the  fixed  ecliptic  is  that  of  1 750,  no  matter  for  what  epoch  the  precession  may 
be  defined.  The  annual  lunisolar  precession  then  diminishes  with  the  time.  But  if 
we  consider  the  precession  upon  the  fixed  ecliptic  of  the  date,  then  the  annual  value 
is  a  quantity  which  increases  with  the  time.  When  we  have  removed  this  ambiguity 
by  adopting  one  or  the  other  of  these  definitions,  we  may  regard  the  lunisolar  pre- 
cession as  dependent  solely  upon  the  motion  of  the  equator,  and  as  defining  that 
motion.  The  same  motion  is  defined  by  the  quantity  w,  which  expresses  its  annual 
value  when  measured  on  a  great  circle  of  the  sphere.  On  the  other  hand,  the  two 
general  precessions  are  resultants  of  both  motions — that  of  the  ecliptic  and  that  of  the 
equator. 

To  avoid  the  confusion  which  may  thus  arise  in  treating  the  subject,  I  have  pro- 
posed the  term  "  Precessional  Constant"  for  a  certain  function  P  of  the  masses  of  the 
Sun,  Earth,  and  Moon,  and  of  the  elements  of  the  orbits  of  the  Earth  and  Moon,  on 
which  the  motion  of  the  equator  depends.  This  quantity  changes  only  through  the 
diminution  of  the  Earth's  eccentricity,  and  the  change  is  so  slight  that  it  can  become 
appreciable  only  after  several  centuries.  From  the  value  of  P  and  the  motion  of  the 
ecliptic,  the  four  motions  of  the  precession  are  derived  by  the  formula?, 

General  precession;  $  —  P  cos  £  —  x  sin  L  cos  e 

Lunisolar  precession ;       pl  r=  P  cos  e 
Precession  in  R.  A. ;         m  —  P  cos2  e  —  x  sin  L  cosec  e 
Precession  in  declination ;  n  —  P  cos  e  sin  £ 

Here  L  is  the  longitude  of  the  instantaneous  axis  of  rotation  of  the  ecliptic,  x  is 
its  annual  motion,  and  £  is  its  obliquity. 

Section  I. 

Previous  determinations. 

Of  the  authoritative  determinations  of  the  precessional  motion  made  during  the 
nineteenth  century  the  first  is  that  of  Bessel.  It  rests  upon  a  comparison  of  the  star 
catalogue  of  the  Fundamenta  Astronomicc  with  the  catalogue  of  Piazzi  for  1800,  the 
latter  being  first  corrected  by  more  recent  Konigsberg  observations.  The  values  of 
the  motions  as  they  are  quoted  in  Engelmann's  edition  of  Bessel's  Abhandlungen 
from  the  Memoirs  of  the  Berlin  Academy  of  Sciences  for  181 5*  are  found  below.  But 

♦BESSEL'S  Abhandlungen,  von  R.  Engelmann,  Vol.  I,  pp.  262-285. 
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in  the  Tabulae  Regiomontanse,  published  in  1830,  somewhat  different  values  are  given. 
The  two  sets  of  values  are  as  follows : 

Abhandlungen.  Tabulae  Regiomontanae. 

//  1/  //  // 

G-en.  prec;  50.18828  +  .02442966T        50.22350 -f  .02443T 
Lunisolar;  50.32832  —  .02435890  50.36354  — .02436 

m  =  46.01 135  +  .03086450         46.04367  +  .03086 
n  —  20.04554  — .00970204  20.05957  — .00970 

T  is  here  counted  from  1 800.0.  in  terms  of  the  century  as  the  unit. 

The  interval  between  the  two  catalogues  on  which  this  determination  rested  was 
only  forty-five  years.  In  1 840,  the  shortness  of  this  period  and  the  fact  that  much  material 
had  accumulated  for  improving  the  values  of  the  constants  in  question,  led  Otto  Struve 
to  make  a  new  determination.  This  determination,  with  some  alterations  by  Peters, 
lias  been  most  in  use  to  the  present  time.  It  depends  upon  a  comparison  of  the 
positions  of  400  stars  in  the  Fundamenfa  with  determinations  made  principally  at 
Dorpat  about  the  mean  epoch  1825.  It  therefore  rests  upon  an  interval  of  about,  seventy 
years*  The  marked  feature  of  this  determination  was  that  it  included  a  determi- 
nation of  the  solar  motion,  which  was  thus  eliminated  from  the  result.  The  corrections 
to  the  provisional  value  of  the  motion  of  precession  for  seventy  years,  which  I  suppose 
to  be  that  of  Bessel,  were  found  to  be: 

In  R.  A. ;  Jm  =  +  1  ".6743 ;  Jp  =  +  1 ".  1 6  ±  o".67 
In  Decl.;  J»=4- o".2624;  Jp  =  +  o".66 ± o".86 
The  combined  result  is  Ap  —  +  o".97  db  ©".53. 
p  —  50". 23492  for  the  epoch  1790. 

In  the  Astronomische  Nachrichten,  No.  485,  Struve  the  elder  quotes  the  result  as 

5o".23449±o".oo77i 

Neither  of  the  Struves  give  any  precessional  motion  but  this.  The  motions  with 
which  Struve's  name  is  associated  are  found  in  the  memoir  of  Peters,  Numerns  Con- 
stans  Nntationis.  Here,  on  page  160,  Peters  gives  50". 3 798  as  the  lunisolar  preces- 
sion for  1800  found  by  Struve.  He  afterwards  deduces  the  following  values  of  the 
various  quantities  on  which  the  precessional  motions  depend,  as  those  resulting  from 
Stbuve's  work: 

//  // 

General  precession,    1800;  50.241 1 -f- .02268T 
Lunisolar  precession,  1 800 ;  50.3 798  —  .02 1 68 

m  —  46.0623  +  .02849 
n  —  20.0607  —  .00863 

These  values  are  tabulated  on  page  200  of  the  paper,  and  are  those  in  most  gen- 
eral use  up  to  the  present  time.    It  will,  however,  be  seen  that  the  value  of  the  general 

♦Memoires  de  l'Academie  Itnpenale  des  Sciences  de  Saint-Pdtersbourg,  sixienie  serie.  Sciences  Matli- 
e'niatiques  et  Physiques,  Tome  III. 
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precession  is  larger  by  o".43  per  century  than  that  found  by  Struve  in  his  printed 
paper.  I  am  not  aware  of  the  origin  of  this  change,  but  its  effects  was  to  increase 
the  error  of  a  determination  already  somewhat  too  large. 

As  I  pointed  out  to  the  Conference  in  1 896,  some  confusion  has  arisen  from  the 
fact  that  in  applying  Peters's  formulae  to  the  trigonometric  reduction  of  star  positions 
the  obliquity  of  the  equator  to  the  fixed  ecliptic  has  been  earned  only  to  terms  in  f, 
whereas,  to  make  the  formula;  consistent  and  homogeneous,  this  quantity  should  be 
carried  to  terms  in  f3. 

Before  the  above  values  had  entirely  superseded  those  of  Bessel,  suspicion  arose 
that  they  were  too  large,  and  that  those  of  Bessel  were  really  nearer  the  truth, 
although  the  data  on  which  they  rested  were  much  inferior.  In  his  "Annales  de 
l'Observatoire,"  Vol.  II,  Le  Verrier  adopts  the  following  values: 

//  // 

General  precession  (1850);  50.23572  +  -022578T 
Lunisolar  precession ;  50.3 7 1 40  —  .02 1 762 

m  —  46.06010  +  .028373 
n—  20.05240  —  .008663 

These  have  been  used  in  the  reductions  of  stars  at  the  Paris  Observatory. 

In  his  "Lehrbuch  eur  Bahnbestimmung  der  Kometeh  und  Plancten,"  Vol.  I  second 
edition,  Oppolzer  develops  the  most  exhaustive  theory  of  precession  which  has  yet 
appeared,  but  without  any  discussion  of  the  observations  on  which  his  adopted  values 

of  the  constants  depend.    His  results  are  (pp.  202,  206) : 

//  // 

General  precession  ( 1 850) ;   50  23465  -f  .022 580T 
Lunisolar  precession ;  50.36924  —  .02 1 7 76 

m  =  46.0593 1  +  .028390 
n  =  20.05 1 5°  —  .008668 

About  1886  Ludwio  Struve  took  advantage  of  the  completion  of  Auwers'  reduc- 
tion of  Bradley's  observations,  and  of  the  Poulkowa  catalogue  for  1845,  1855,  and 
1865,  to  enter  upon  the  most  thorough  discussion  of  the  observations  which  has  yet 
been  published.  The  negative  correction  found  to  Otto  Struve's  result  was  quite 
large.    The  precessional  motions  m  and  n,  found  by  L.  Struve,  were 

//  // 

General  precession  (1850);  50.2283  +  .02120T 
Lunisolar  precession ;  50.3407  —  .02132 

m  =  46.05  54  +  .02  74 1 
n  —  20.0452  —  .00849 

It  must  be  stated,  however,  that  the  smallness  of  these  values  arises  partly  from 
the  large  positive  correction  which  Auwers  applied  to  Bradley's  Right  Ascensions. 
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When  the  various  determinations  are  reduced  to  a  common  equinox,  it  is  found  that 
the  discordances  are  materially  diminished.* 

All  the  preceding  determinations  rest  upon  Bradley  stars,  which  I  conceive  to 
afford  the  best  material  available  at  the  present  time.  Determinations  from  a  great 
number  of  faint  stars  have  been  made  by  Dreyer,  Nyren,  and  Bolte. 

Dbeyer's  investigation  is  found  in  the  Journal  Copernicus,  Volume  II,  Dublin, 
1882.  It  is  prefaced  by  a  very  exhaustive  history  of  previous  determinations  of  this 
constant.  The  material  consists  of  Lalande's  catalogue,  and  the  catalogue  of  Schjel- 
lerup  for  1865.  The  system  of  right  ascensions  to  which  both  catalogues  are  reduced 
is  N,.   The  results  for  1800  + 1  are— 

//  // 

Lunisolar  precession  —  50. 3752  —  0.0002 1 68/ 
General  precession  =50.2365+0.0002268/ 
m  —  46.058 1  +  0.0002849/ 
n  —  20.0589  —  0.0000862/ 

Nyren's  result  was  based  on  a  comparison  of  stars  in  Bessel  zones  with  Schjel- 
lerup,  and  is  much  smaller  than  the  above.  The  motion  is,  however,  greatly  increased 
when  the  two  catalogues  are  reduced  to  the  system  N(,  so  that,  ultimately,  the  discord- 
ance is  not  great.    But  the  brevity  of  the  interval  renders  the  determination  weak. 

Bolte's  investigation  is  based  on  the  declinations  of  Lalande  as  compared  with 
Schjellerup.  His  work  was  published  only  as  a  doctorial  thesis,  and  is  therefore  not 
so  well  known  as  it  should  be. 

These  determinations  would  be  better  than  any  that  could  be  made  from  Bradley 
stars  were  it  not  for  the  probable  systematic  errors  in  the  observations  and  reductions. 
First  among  these  errors  is  that  in  Right  Ascension  depending  on  the  magnitude  of 
the  star.  When  we  come  down  to  faint  stars  this  error  may  increase  to  an  almost 
unknown  extent  Moreover,  as  a  general  remark  applicable  to  all  zone  observations, 
it  may  be  said  that  there  is  great  room  for  systematic  errors,  both  in  the  Right 
Ascension  and  Declination,  which  would  materially  affect  the  value  of  the  precessional 
motion  derived  from  them. 

List  of  the  principal  determinations. 

For  the  sake  of  easy  comparison  the  values  above  given,  with  some  others,  are 
reduced  to  1850,  and  the  results  shown  below. 

In  the  reductions  of  the  earlier  values  to  1 850  the  motion  in  my  "  Elements  and 
Constants  "  has  been  used. 

The  value  N0  is  the  preliminary  one,  used  in  the  tables  of  the  planets  and  in  my 
"  Elements  and  Constants." 

*See  Vierteljahrsschrift  der  Astronomischen  Gesellschaft,  Vol.  XIII,  187S,  page  107. 
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The  value  N,g  is  that  reached  in  the  present  paper. 
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Section  II. 

Material  used  in  the  present  discussion. 

In  starting  with  a  new  determination,  the  first  question  to  be  considered  is  what 
material  we  shall  use.  What  we  want  is  the  annual  motion  of  the  equinox,  arising 
from  the  combined  motions  of  the  equator  and  the  ecliptic,  relative  to  an  axis  abso- 
lutely fixed  in  space.  As  observations  can  not  be  referred  to  any  line  or  plane  which 
we  know  to  be  absolutely  fixed,  we  are  obliged  to  assume  the  general  mean  position 
of  the  fixed  stars  as  affording  the  invariable  axes  of  reference.  Here  we  meet  with 
difficulties  arising  from  the  proper  motions  of  the  stars,  and  especially  from  the  solar 
motion.  The  elimination  of  the  latter  from  the  mean  result  is  a  matter  of  especial 
difficulty. 

It  is  also  to  be  recalled  that  the  position  of  the  equinox  relative  to  the  mean 
stellar  sphere,  assumed  as  fixed,  requires  three  practically  separate  determinations. 
These  are — 

1.  The  correction  to  the  Sun's  absolute  mean  right  ascension  or  longitude.  This 
is  obtained  principally  from  observations  of  the  Sun's  declination. 

2.  The  correction  to  the  general  mean  right  ascension  of  the  clock  or  standard 
stars  relative  to  the  Sun. 

3.  The  determination  of  the  positions  of  the  clock  stars  relative  to  the  great  mass 
of  stars  assumed  to  define  the  fixed  celestial  sphere. 

In  the  present  investigation  it  is  needless  to  consider  the  first  two  determinations 
because,  by  a  resolution  of  the  Conference  of  1 896,  in  accordance  with  which  this 
investigation  is  made,  the  position  of  the  equinox  among  the  clock  stars  is  defined  by 
the  system  N„  formed  by  those  equatorial  stars  whose  position  is  given  in  the  Appen- 
dix to  the  Catalogue  of  1098  Clack  and  Zodiacal  Stars,  found  in  Volume  I  of  these 
Astronomical  Papers.  It  is  therefore  only  necessary,  so  far  as  the  Right  Ascensions 
are  concerned,  to  determine  the  positions  of  the  great  mass  of  stars  relative  to  the 
standard  stars  in  question. 
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The  question  now  arises,  what  stars  we  shall  choose.  Undoubtedly,  the  great 
mass  of  the  fainter  stars,  down  to  the  ninth  magnitude,  which  we  find  in  various  zones 
observed  during  the  present  century,  would  be  preferable.  But  a  determination  made 
from  these  stars  would  be  subject  to  very  large  systematic  error,  arising  from  the 
shortness  of  the  interval  between  the  earlier  and  the  later  observations.  This  uncer- 
tainty is  magnified  by  the  possibility  of  very  large  personal  error  due  to  the  small 
magnitude  of  the  star.  It  is  now  well  known  that,  at  the  present  time,  especially 
when  observations  are  made  by  the  chronograph,  faint  stars  are  observed  too  late  rela- 
tively to  brighter  ones.  Were  this  personal  error  constant  at  all  times  and  for  all 
observers  during  the  century,  no  error  would  thus  arise.  But  there  is  strong  reason 
for  suspecting  that  this  error  not  only  varies  from  one  observer  to  another,  but  that  it 
is  much  larger  for  the  recent  than  for  the  older  observations.  If  such  is  the  case,  the 
precessional  motion  derived  from  such  a  comparison  would  be  too  great.  The  proba- 
bility of  this  error  has  led  me  to  depend  mainly  upon  observations  of  the  brighter 
stars,  especially  those  in  Bradley's  catalogue.  As  the  total  number  of  stars  exceeds 
3,000,  it  may  be  expected  that  all  accidental  errors  will  be  nearly  eliminated  from  the 
mean.  Though  the  systematic  error  just  considered  may  not  be  wholly  eliminated,  it 
will  certainly  be  much  smaller  than  in  the  case  of  a  determination  founded  on  much 
fainter  stars  observed  during  a  shorter  interval. 

I  have,  however,  made  an  exception  of  the  Lalande  catalogue  by  using  the  com- 
parisons of  this  catalogue  with  those  of  Boss  and  Schjellerup.  The  former  deter- 
mined his  personal  error  dependent  on  magnitude.  In  the  case  of  Schjelleri/p  this 
error  seems  to  be  small. 

Section  III. 

Provisional  2^eccssions. 

From  the  considerations  already  set  forth,  the  desirableness  will  be  seen  of  adopt- 
ing at  the  start  a  provisional  value  of  the  precessional  constant  which  will  be  as  near 
as  practicable  to  the  truth.  That  Struve's  value  needs  a  large  negative  correction  is 
now  widely  admitted,  and,  it  seems  to  the  writer,  made  quite  clear  by  the  data  already 
set  forth  in  the  preceding  pages.  I  shall,  therefore,  take  as  provisional  values  the 
precessional  motion  already  alluded  to,  which  is  employed  in  my  "Elements  and  Con- 
stants," as  the  provisional  value  to  be  corrected. 

The  principal  quantities  for  the  epoch  1850  are  these: 

The  precessional  constant  for  a  solar  century  of  36524.2  days 

P  =  5489".  78  -  o".oo364T 

// 

General  precession;  p  —  5023.71 
Lunisolar  nreeession ;  pt  —  5036.02 

100  mi  —  4606.36 

100  n  —  2004.79 

The  fundamental  quantity  to  be  determined  is  the  correction  to  this  value  of  P; 
but  it  will  be  more  convenient  to  take  pt,  the  lunisolar  precession,  sis  the  Quantity  to  be 


12 


THE  PRECESSIONAL  CONSTANT. 


corrected.  The  correction  I  shall  call  Sp,  or  simple  p,  so  that  the  true  lunisolar  pre- 
cession sought  shall  be 

p,  =  5036"  02  +  SP,  or  5°36"-°2  +P 
Section  IV. 
Systematic  corrections  to  the  proper  motions. 

I  have  taken,  as  the  basis  of  the  work,  Auwebs'  proper  motions  of  the  Bradley 
Stars,  as  found  in  Vol.  Ill,  Abth.  II  of  his  reduction  of  Bradley's  observations.  But 
it  is  necessary  to  correct  the  proper  motions  of  Auwers  so  as  to  reduce  them  to  the 
system  already  employed  for  the  fundamental  stars.  These  motions,  it  will  be  recalled, 
are  derived  by  a  comparison  of  the  Bradley  catalogue  for  1755  with  a  number  of  the 
best  recent  determinations.  The  latter  were  reduced  to  a  mean  of  the  Greenwich 
catalogues  for  i860  and  1864.  The  difference  between  the  mean  corrections  of  these 
catalogues  and  of  the  Bradley  catalogue  will  therefore  be  the  correction  to  the  proper 
motions  for  about  one  hundred  and  five  years. 

In  the  Declinations,  however,  we  have  a  simpler  way  of  proceeding.  The  proper 
motions  in  question  were  incorporated  in  the  A.  G.  C.  fundamental  catalogue,  Publication 
XIV  of  the  Astronomische  Gesellschaft,  and  may  therefore  be  considered  as  system- 
atically identical  with  these.  Now,  in  the  work  on  the  fundamental  stars,  I  have 
determined  the  correction  to  reduce  the  declinations  of  the  A.  G.  C.  to  those  of  the 
normal  system,  together  with  their  centennial  variations.  The  latter  are  clearly  the  cor- 
rections to  be  applied  to  the  proper  motions  of  Auwers-Bradley  in  Declination. 

The  systematic  corrections  to  the  proper  motions  in  question  thus  consist  of  two 
classes;  those  for  reduction  to  the  normal  system,  and  those  for  reduction  to  the  new 
provisional  value  N0  of  the  precessional  constant. 

It  is  found  that,  for  reduction  to  the  system  N, ,  the  Right  Ascensions  of  Bradley 
required  the  constant  correction  for  equinox 

Ja  1755  —  —  0*079 

The  modern  standards  to  which  the  Right  Ascensions  are  reduced  in  publication 
XIV,  and  in  Auwers-Bradley,  is  that  of  the  two  Greenwich  seven-year  catalogues 
for  1 860  and  1 864,  the  Right  Ascensions  of  the  second  one,  that  for  1 864,  being  first 
increased  by  o'.oi.  By  the  comparison  of  these  catalogues  with  the  system  N,  it 
appears  that,  this  system  would  require  the  following  constant  correction  for  equinox: 

Catalogue  for  1 860  N,  —  A  =  N,  —  Grh.  =  +  o".024 

Catalogue  for  1 864  N,  —  A  —  N,  —  Grh.  —  o'.oi  =  -f  o".024 

There  is,  therefore,  an  exact  agreement  between  the  results  derived  by  comparison 
with  the  two  catalogues.  Thus  the  constant  corrections  to  the  A.  G.  system  on  account 
of  equinox  are 

s. 

Epoch  1755  —-079 

Epoch  i860  +  024 

Correction  for  105  years,  .  .  +.103 
Correction  for  100  years,    .    .  +.098 

a  result  which  agrees  exactly  with  that  found  by  direct  comparison  of  N,  with  A.  G.  C 
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We  must  also  apply  the  mean  periodic  term  for  the  two  Greenwich  Catalogues,  of 
which  the  value  is 

da  —  +  *.oi$  cos  a  —  •.oil  sin  ar='.oi8  cos  (a  +  380) 

If  to  the  constant  '.098  we  apply  the  difference  of  the  corrections  Jat  for  the  two 
catalogues,  we  shall  liave  the  following  values  of  the  systematic  corrections  to 
Auwebs'  proper  motions  in  Right  Ascension  in  order  to  reduce  them  to  the  modern 
system  and  to  Struve's  constant. 


Dec. 

Corr. 

0 

s. 

-30 

+  .140 

—  20 

.129 

— 10 

.126 

—  0 

.114 

+  10 

.097 

20 

.082 

30 

•075 

40 

.102 

50 

.096 

60 

.14 

70 

.16 

80 

.1.6 

90 

.16 

For  the  precessions  at  the  mid-epoch  18 10,  Auweks  uses 

//  // 

m  —  4606.5 1 ;  n  —  2005.98 
N0  has  4605.24  2005.13 
A-N0=  +1.27  +0.85 

Hence,  the  centennial  proper  motions  of  Auwers  in  Right  Ascension  require  the  cor- 
rection 

//  // 
+  1.27  -I-0.85    sin  a  tan  6 
s.  s. 

=  0.085+0.057  sin  a  tan  6 

In  the  equations  of  condition  the  proper  motions  are  multiplied  by  cosine  8.  If 
we  put  C  for  the  numbers  given  above  as  a  function  of  the  Declination  increased  by 
i".27,  the  complete  reduction  required  is  given  by  the  expression 

da  cos  <J  =  C  cos  <J  +  o".27  cos  8  cos  (a  +  380)  +  o".85  sin  a  sin  8 

This  quantity  has  been  tabulated  in  a  table  of  double  entry  as  a  function  of  both 
coordinates. 

In  the  case  of  the  Declinations  the  corresponding  corrections  to  the  proper 
motions  will  consist  of  the  following  three  parts : 

1 .  The  coefficient  of  T  in  the  reduction  of  the  A.  G.  to  the  new  normal  system 
CM,  which  correction  is  a  function  of  the  Declination  alone. 
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2.  The  negative  of  the  correction  JSa  to  Auwers-Buadley,  which  will  be  tabulated 
subsequently  in  Part  II  as  a  function  of  the  Right  Ascension  and  Declination. 

3.  The  reduction  for  the  precession  N„,  of  which  the  value  is  4^'  —  +  o".85  cos  nr. 
The  sum  of  these  three  quantities  also  was  tabulated  in  a  table  of  double  entry. 
The  tables  of  the  preceding  corrections  constructed  for  use  in  the  work  are  not 

printed  in  the  present  paper. 

Section  V. 

Formation  and  treatment  of  tlie  equations  of  condition. 

The  fundamental  quantity  of  which  the  correction  is  to  be  found  is  that  already 
defined  as  the  precessional  constant.  But  we  may  equally  take  for  correction  the 
value  of  the  lunisolar  precession 

Px  =  P  cos  e 

as  the  quantity  for  which  the  correction  is  to  be  determined,  because  the  value  of  the 
obliquity  £  is  so  well  known  that  the  correction  of  pu  may  be  regarded  as  a  constant. 
We  shall  then  have 

m  —  jh  cos  e  —  x  sin  L  cosec  e 
n=pi  sin  e 

The  planetary  precession,  u  sin  L  cosec  e,  can  be  derived  only  from  theory,  and  the 
value  already  found  from  the  masses  of  the  planets  must  be  adhered  to.  We  there- 
fore have  for  the  symbolic  corrections  to  m  and  n 

8  m  —  cos  e  8  p 
8  n  —  sin  e  8  p 

The  annual  precession  in  R.  A.  being 

m  +  n  sin  a  tan  8 

its  symbolic  correction  is 

8  p  (cos  £  -+- sm  e  8m  a  t*111  *0 

If  we  put  fx  and  //'  for  the  proper  motions  in  Right  Ascension  and  Declination  as 
derived  with  the  provisional  value  of  the  precession,  the  equations  of  correction  will 
be  those  which  will  reduce  this  proper  motion  to  zero.  Omitting  for  the  moment  the 
solar  motion,  these  equations  will  be 

8  p  (cos  £  cos  8  +  sin  £  sin  a  sin  8)  —  fx  cos  8 
8  p  sin  £  cos  a  =  n' 

If  we  add  to  our  equations  the  terms  depending  on  the  coordinates  X,  Y,  Z  of 
solar  motion,  the  complete  equations  of  condition  will  be 

,  sin  a  v     cos  a  _r  (1) 
aSP  +  ~~y  X  p~  Y  =  P  008  6 

,  „  cos  a  sin  8  sin  a  sin  8  cos  8  , 
°&P  +  2  X  +  Y-  —  Z  =  M' 
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where  we  have  put 

a  =  cos  e  cos  8  -f-  sin  e  sin  8  sin  a 
b  rr  sin  c  cos  a 

These  equations  are  well-known  ones  heretofore  used  by  others.  I  have  omitted 
all  expressions  representing'  a  possible  revolution  of  the  stellar  system,  believing  that 
their  introduction  would  be  useless. 

We  shall  have  as  many  pairs  of  equations  of  this  kind  as  there  are  stars  whose 
proper  motions  are  used.  From  these  equations  we  are  to  derive  the  values  of  the  tliree 
components  of  solar  motion  X,  Y,  and  Z,  and  the  correction  of  the  lunisolar  preces- 
sion 8 p.  But  in  proceeding  to  subject  these  equations  to  the  method  of  least  squares, 
we  meet  with  two  serious  and  embarrassing  difficulties.  The  first  of  these  arises  from 
the  presence  in  each  equation  of  the  quantity  p,  the  distance  of  the  star.  We  have 
no  criterion  for  the  value  of  p  except  the  proper  motion  of  the  star  and  its  magnitude. 
Both  of  these  criteria  are  extremely  uncertain,  owing  to  the  great  diversity  in  the 
actual  magnitudes  of  the  different  stars  and  the  actual  amount  of  their  proper  motions. 

A  very  serious  objection  to  using  the  proper  motion  of  a  star  as  a  criterion  of  its 
distance  is  that  we  thereby  assume  the  solution  of  the  very  problem  at  which  we  are 
aiming.  The  amount  of  the  computed  proper  motion  is  a  function  of  the  assumed 
value  of  the  processional  constant,  and  is  also  affected  by  the  solar  motion,  which  is 
inseparably  connected  with  the  determination  of  the  precessional  constant.  If,  then, 
we  select  those  stars  which  have  very  small  apparent  proper  motions,  our  selection 
will  not  consist  of  those  whose  motions  are  absolutely  the  smallest,  but  of  those  which 
give  small  corrections  to  the  precessional  constant  and  to  the  adopted  value  of  the 
solar  motion.  Thus  will  arise  a  bias  in  the  direction  of  the  values  with  which  we 
start.  I  therefore  conclude  that,  while  it  will  be  advisable  ultimately  to  assign 
weights  dependent  011  the  amount  of  the  proper  motion,  we  can  do  this  only  with  the 
residual  proper  motions  which  remain  after  making  the  best  allowance  we  can  for  the 
solar  motion  and  for  the  correction  to  the  precessional  constant. 

This  reasoning  applies  mainly  to  the  cases  of  individual  stars.  If  we  consider 
large  classes  of  stars,  the  force  of  the  objection  will  be  greatly  diminished,  since  the 
mean  amount  of  the  proper  motion  will  probably  then  approximate  to  the  inverse  of 
the  mean  distance.  We  shall  therefore  follow  previous  investigators  in  classifying  the 
stars  according  to  their  magnitudes.  We  shall  then  estimate  the  mean  distance,  or, 
rather,  the  mean  parallax,  of  the  stars  of  each  class  to  be  proportional  to  the  mean 
value  of  the  proper  motion  for  the  individual  stars  of  the  class.  This  proceeding  will 
have  the  advantage  of  enabling  us  to  make  a  combination  by  which  the  uncertainty 
as  to  the  mean  value  of  p  in  each  case  will  be  merged  with  the  question  of  the  acci- 
dental errors,  as  follows: 

Let  us  multiply  the  preceding  equations  by  p.  The  equations  of  condition  (1) 
will  then  take  the  following  form: 

a  p  8  p  +  sin  or  X  —  cos  aY  —  /*  P  cos  8  —  pn 

bp  8 p-\-  cos  a  sin  8  X  +  sin  a  sin  8  Y  —  cos  8  Z—  fx'  p  —  pn'  ' 

It  is  quite  indifferent  in  writing  the  equations  whether  we  consider  the  latter  in 
this  or  their  original  form,  since  we  shall  have  in  either  case  to  multiply  each  equation 
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by  such  factor  or  weight  as  will  reduce  the  probable  errors  of  the  second  term  to  a 
constant  value  for  all  the  equations.  If  in  using  the  first  or  original  form  we  have 
the  system  of  weights 

u\;  w3;  etc., 

it  is  evident  that  in  using  the  second  form  our  system  of  weights  should  be 

Pi'      P*  P* 

We  shall  then  have  the  same  equation,  whether  we  start  with  one  form  or  the  other 
In  either  case  the  weight  to  be  chosen  will  be  that  which  will  reduce  the  probable 
eiTor  of  the  absolute  terms  to  the  same  value  in  the  case  of  each  group  of  stars. 
This  probable  error  is  composed  of  two  parts — the  absolute  proper  motions  of  the  stars 
themselves,  as  seen  from  the  earth,  and  the  errors  of  their  determinations.  But  it  is 
unnecessary  to  consider  these  two  sources  of  error  separately,  because  their  combined 
effect  is  shown  by  the  mean  divergence  of  the  individual  stars  from  the  general  mean 
in  the  case  of  each  class  of  stars. 

It  follows  that  if  we  make  the  proposed  classification,  and  find  the  mean  value  of 
the  residual  proper  motion  in  the  case  of  the  stars  of  each  class,  the  square  roots 
of  these  means  will  form  a  series  of  divisors  for  the  several  equations  of  condition 
which  will  reduce  the  probable  errors  of  all  to  the  same  value. 

This  difficulty  disposed  of,  there  remains  one  yet  more  embarrassing,  arising  from 
the  fact  that  the  proper  motions,  which  we  are  compelled  to  treat  as  residual  errors, 
do  not  follow  the  normal  law  of  error,  even  when  reduced  to  one  standard  of  mean 
parallax  in  the  way  proposed 

In  all  determinations  of  the  processional  motion  heretofore  made,  the  proper 
motions  of  the  stars  have  been  treated  as  residual  errors  by  the  application  of  the 
method  of  least  squares.  The  acceptance  of  a  mean  of  several  different  observed 
values  of  the  same  quantity  as  the  most  likely  value  is  a  special  example  of  this  method. 
But  the  method  has  to  be  modified,  in  so  far  that  stars  of  excessively  large  proper 
motion  are  arbitrarily  thrown  out.  Ludwig  Struve,  in  his  discussion  on  the  subject, 
adopted  the  principle  of  throwing  out  those  stars  whose  proper  motion  was  more  than 
ten  times  the  mean  proper  motion  of  stars  having  that  magnitude.  But  it  seems  that 
this  led  to  the  rejection  of  only  seven  stars  in  all. 

That  this  system  of  rejecting  all  results  which  diverge  largely  from  the  general 
mean  is  logically  sound  has  long  been  understood.  But  there  is  an  evident  want  of 
system  in  the  proceeding,  because  we  assume  a  certain  fixed  limit  within  which  results 
are  accepted  with  their  full  weight,  while  without  the  limit  they  are  totally  rejected. 
The  quantities  which  deviate  most  largely  are  those  whose  rejection  or  retention  will 
most  influence  the  concluded  result.  Consequently,  the  final  result  will  frequently 
depend  largely  on  where  we  set  the  limit  of  rejection,  and  this  is  a  matter  which  has 
to  be  left  entirely  to  the  judgment  of  the  investigator.  I  have  on  several  occasions 
pointed  out  that  the  cases  in  which  the  m«thod  of  least  squares  gives  the  best  result  are 
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not  the  typical  ones  of  astronomy.*  The  general  propositions  which  relate  to  the 
subject  may  be  formulated  as  follows: 

I.  The  method  of  least  squares  gives  the  best  value  when  and  only  when  the 
probability  of  making  an  error  of  a  given  amount  x  is  expressible  as  a  function  of  x 
by  the  formulae, 

II.  In  astronomical  practice  it  commoniy  happens  that  the  actual  probability  of 
an  error  of  given  magnitude  x  is  not  given  by  a  function  of  this  form,  but  deviates 
from  this  function  in  such  a  way  that,  when  the  closest  adjustment  to  the  formula  is 
made,  small  errors  and  large  errors  are  more  frequent  than  intermediate  ones.  • 

III.  In  this  case  the  best  result  is  obtained  by  giving  to  the  more  divergent  results 
a  weight  depending  on  the  magnitude  of  the  divergence,  gradually  converging  toward 
zero  as  the  divergence  increases. 

IV.  The  necessity  of  thus  modifying  the  method  of  least  squares  becomes  greater 
the  wider  the  deviation  of  the  actual  distribution  of  error  from  that  given  by  the 
exponential  formulas. 

In  order  to  judge  of  the  best  mode  of  procedure,  I  begin  by  arranging  the 
reduced  proper  motions  of  the  Bradley  stars  in  the  direction  of  R.  A. — that  is,  the 
values  of  /<  cos  6  derived  by  the  methods  already  described,  in  such  a  way  as  to 
show  their  distribution  statistically.  In  order  to  separate  the  effect  of  the  parallactic 
motion,  I  have  made  the  classification  not  only  according  to  magnitude,  but  according 
to  the  quadrant  of  right  ascension.  The  adopted  magnitudes  are  those  given  by 
Auwers  in  the  catalogue.  As  the  stars  of  the  first  magnitude  are  too  few  to  be  used 
by  themselves,  I  have  put  all  down  to  magnitude  2.9  in  one  group.  A  second  group 
includes  those  from  3.0  to  3.9 ;  a  third  those  of  each  following  magnitude  down  to 
6.9.  All  below  this  magnitude  I  have  included  in  one  group  as  of  magnitude  7.  For 
each  quadrant  of  right  ascension  and  each  group  I  have  taken  the  proper  motions  in 
right  ascension  and  made  a  count  of  the  number  having  certain  values.  One  of  the 
subgroups  thus  formed  comprises  the  proper  motions  which  come  out  as  exactly  zero. 
The  negative  proper  motions  are  then  classified  as  —  o".  and  a  fraction,  —  1".  and  a 
fraction,  etc.,  up  to  —  16".  and  a  fraction.  All  of  if'.o  and  higher  are  put  into  one 
group.  A  similar  division  is  made  on  the  positive  side.  The  result  of  this  process  is 
shown  in  Table  I. 

Whatever  law  of  probability  of  error  we  adopt,  there  are  always  two  quantities 
expressive  of  the  mean  magnitude  of  the  error  wliich  admit  of  rigorous  definition. 
One  of  these  is  the  mean,  the  other  the  probable  error.  The  mean  of  the  squares  of 
all  the  errors  made  is  evidently  a  quantity  which  can  be  determined  in  all  cases. 
Moreover,  by  actual  enumeration  we  can  always  have  two  limits  between  which  one- 
half  the  actual  deviations  fall,  and  we  can  choose  those  limits  so  that  their  difference 
shall  be  a  minimum.  We  may  regard  one-half  of  this  difference  as  the  probable 
error.    We  are  not,  however,  confined  to  one-half  as  a  fraction.    We  may  take  any 

'American  Journal  of  Mathematics,  Volume  VIII,  page  343,  and  Elements  and  Constants,  Chapter  IV. 
A  P,  VOL  VIII,  PT  I  a 
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other  fraction  of  the  whole  number,  and  find  the  closest  limits  between  which  that 
fraction  is  contained.  The  difference  of  these  limits  will  in  all  cases  be  a  fairly  deter- 
minate quantity. 

Table  I. — Distribution  in  magnitude  of  the  proper  motions  n  cos  6  of  the  Bradley  Stars  in  R.  A. 
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Section  VI. 

Statistical  method  of  determining  the  precession. 

I  have  used  the  numbers  in  Table  I  to  ascertain  what  value  of  the  precession 
in  right  ascension  would  result  when  we  do  not  take  means,  but  simply  assume  that 
negative  and  positive  proper  motions  are  equally  probable,  and  that  the  most  likely 
value  of  the  precession  is  that  which  makes  the  number  of  positive  and  of  negative 
motions  equal.  This  determination  can  be  made  by  actual  count  for  each  magnitude. 
But  if  we  proceeded  thus  without  respect  to  the  quadrant  of  right  ascension,  we  should 
arrive  at  an  incorrect  result  on  account  of  the  effect  of  solar  motion,  combined  with 
the  unequal  distribution  of  the  stars.    For  example,  if  there  are  an  excess  of  stars  in 
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that  semicircle  of  right  ascension  in  which  the  effects  of  the  solar  motion  are  positive, 
it  is  evident  that  the  result  would  thus  be  biased.  We  shall  therefore  determine,  for 
each  quadrant  of  right  ascension,  what  value  of  the  proper  motion  will  correspond  to 
that  quadrant  in  such  away  that  the  number  of  positive  and  negative  excesses  of  actual 
proper  motion  shall  be  equal.  This  value  I  call  the  mid-point.  The  method  of  find- 
ing it  in  the  case  of  each  quadrant  and  each  magnitude  is  shown  in  the  following 
tables.  In  any  one  case  is  given  first  the  number  of  negative  proper  motions  and  then 
the  number  of  positive  ones.  We  next  determine  how  far  we  must  count  in  the  direc- 
tion of  the  excess  in  order  to  find  a  number  of  proper  motions  equal  to  one-half  of 
that  excess.  In  finding  this  point  the  proper  motions  within  each  space  of  one  second 
have  been  assumed  as  equally  distributed  throughout  that  second. 

Thus  we  find  the  mid-point  for  each  quadrant  and  each  magnitude  as  shown  below: 

Magnitude  1-2.9  Magnitude  3.0-3.9 

No.  of  proper  motions.  No.  of  proper  motions. 


Quadr. 

Negative, 

Positive. 

Mid-point. 
// 

Negative. 

Positive. 

Mid-point. 
// 

1 

3 

21 

+  3-o 

12 

38 

+  2.7 

TT 
11 

10 

7 

— 1.0 

29 

IO 

5-1 

TTT 
111 

20 

I  I 

-2.5 

37 

19 

—  3° 

TV 
1 V 

2 

12 

+  6.5 

1 1 

31 

-r  34 

Sum 

35 

51 

+  6.0 

89 

98 

—  2.0 

Magnitude  4.0-4.9 

Magnitude  5.0-5.9 

No.  of  proper  motions. 

No.  of  proper  motions. 

Quadr. 

Negative. 

Positive. 

Mid-point. 
// 

Negative. 

Positive. 

Mid-point. 
ft 

I 

35 

8l 

+  i-93 

70 

179 

+  2.O5 

II 

52 

28 

- 1-45 

Il8 

76 

—  O.9O 

III 

55 

35 

—  0.87 

128 

62 

—  2.20 

IV 

36 

105 

+  2.1 1 

66 

234 

+  2.25 

Sum 

178 

249 

+  1.72 

382 

551 

+  I.20 

Magnitude  6.0-6.9 

Magnitude  7.0  -f 

No.  of  proper  motions. 

No.  of  proper  motions. 

Quadr. 

Negative. 

Positive. 

Mid-point. 
// 

Negative. 

Positive. 

Mid-point. 
tt 

I 

93 

271 

+  2.00 

23 

6l 

+  2.8 

II 

198 

98 

-1-25 

56 

27 

—  1.4 

III 

158 

80 

—  1.94 

18 

16 

—  0.3 

IV 

103 

238 

+  I.70 

23 

51 

+  1.6 

Sum 

552 

687 

+  0.5I 

I20 

155 

+  2.7 

The  mean  of  the  four  mid-points  of  any  one  group  is  the  mid-point  as  it  would 
be  found  were  there  no  parallactic  motion,  and  expresses  the  correction  to  the  mean 
precession  of  all  the  stars  in  R.  A.    The  result  in  the  case  of  each  magnitude  for  the 
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correction  to  the  precession  in  right  ascension  for  a  century,  together  with  the  combi- 
nation of  the  six  results,  is  as  follows : 

Mag.  1-2 
3 
4 
5 
6 

7 


wt. 

+  i-5 

0.1 

—  °-5 

°-5 

+  043 

1.8 

+  0.30 

6-3 

+  0.13 

12.0 

+  0.68 

3-5 

Mean:  +0.27 

I  have  combined  the  results  of  the  several  groups  in  yet  another  way  by  multi- 
plying the  preceding  numbers  of  proper  motions  by  0.5  for  magnitude  1  —  2.9 ;  0.7  for 
magnitude  3 ;  0.9  for  magnitude  4 ;  and  unity  for  the  remaining  three  magnitudes. 
The  result  is  as  follows : 


Quad. 

Negative. 

Positive. 

Mid-point 
// 

I 

227 

621 

+  2.1 1 

II 

444 

237 

-1-25 

III 

390 

208 

-1-73 

IV 

233 

646 

+  2.08 

1294 

1712 

+  1. 21 

Mean  corr.  to  cos  6Dta  —  +  o" .30 

The  two  methods  are  in  fair  accordance,  and  lead  to  the  conclusion  that  by  the 
method  in  question  the  mean  excess  of  the  observed  motions  in  the  direction  of  R.  A. 
above  the  precession  N0  is 

+  o".285 

This  mean  excess  of  the  motion  in  R.  A.  is  really  the  mean  value  of  dm  cos  8. 
It  must  therefore  be  divided  by  the  mean  value  of  cos  6  for  that  part  of  the  celestial 
sphere  covered  by  the  Bradley  stars,  say  from  the  pole  to  20 0  S.  Dec.  This  mean 
we  find  to  be  by  integration 

[cos  5]  =  0.84 

We  thus  have 

A   _o".285      .  „ 

and 

dp  —  dm  sec  «  =  +  o"-37 

as  the  result  of  the  statistical  method. 

The  same  numbers  will  enable  us  to  determine  the  R  'A.  of  the  apex  of  the  solar 
motion.    If  we  put 

q;  the  parallactic  motion  of  a  star  900  from  the  apex  of  the  solar  way; 
A,  D ;  the  R.  A.  and  Dec.  of  that  apex ; 
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In  the  four  lines  of  each  part  of  the  table  for  each  quadrant  are  given  the  numbers 
of  the  motions  for  each  of  the  four  magnitudes,  4,  5,  6,  and  7,  together  with  their  sum. 


Table  II. — Summation  of  distributed  proper  motions  jj  cos  6  for  Mags.  4  to  7. 
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The  influence  of  the  parallactic  motion  is  shown  very  strikingly  by  the  diversity 
of  the  extreme  numbers  in  the  four  quadrants.  It  should  evidently  be  allowed  for  in 
seeking  the  law  of  distribution.    I  therefore  proceeded  as  follows: 

Owing  to  the  narrow  limits  between  which  the  proper  motions  of  each  column  are 
contained,  we  may  regard  the  number  of  proper  motions  given  in  each  column  as 
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corresponding  to  a  definite  proper  motion  numerically  greater  by  o".5  than  the  num- 
ber given  at  the  top  of  the  column.  Bearing  this  in  mind,  we  seek,  by  inspection  of 
the  sums  in  any  one  quadrant,  and  in  each  column,  to  find  the  apex  of  the  curve 
whose  abscissa  are  the  proper  motions  and  whose  ordinates  are  the  number  of  stars 
having  that  proper  motion.  I  have  not  attempted  to  locate  this  apex  with  any  greater 
precision  than  as  being  either  in  each  column  or  midway  between  two  adjacent  col- 
umns. In  the  location  I  have  been  governed  very  largely  by  the  statistical  mid-points 
already  found  for  the  several  quadrants.  Thus  I  have  been  led  to  take — 
In  Quadrant  I,  apex  between  columns  1".  and  2".,  corresponding  to  fx  cos  6  —  +  2".o; 
In  Quadrant  II,  apex  between  columns  —  o".  and  —  1".,  corresp'g  to  ju  cos  d  =  —  i".o; 
In  Quadrant  III,  apex  in  column  —  1".,  corresponding  to  fi  cos  6  =  —  i".5; 
In  Quadrant  IV,  apex  in  column  +  1".,  corresponding  to  fi  cos  6  —  +  i".5. 

From  the  data  thus  found  I  have  formed  the  table  below,  in  Avhich  the  headings 
of  the  columns  are  no  longer  the  proper  motions  themselves  but,  in  the  case  of  each 
quadrant,  the  proper  motion  diminished  by  the  value  corresponding  to  the  apex  of  the 
curve  just  given.  When,  as  in  Quadrants  III  and  IV,  the  apex  is  in  one  of  the  col- 
umns, I  have  taken  the  numbers  as  they  stand  in  Table  II,  as  will  readily  be  seen  by 
inspection.  In  the  column  of  Quadrants  I  and  II,  I  have  taken  the  half  sum  of  the 
corresponding  numbers  in  each  pair  of  adjoining  columns.  Moreover,  I  have  reversed 
the  arrangement  of  the  numbers  so  that  positive  proper  motions  at  the  top  of  the  col- 
umn always  correspond  to  the  direction  of  the  parallactic  motion,  whether  this  is  posi- 
tive or  negative.  Since  this  motion  is  directed  very  nearly  from  R.  A.  =  2700,  the 
numbers  of  the  second  and  third  quadrants  are  reversed  in  adding  them  to  the  others. 

Table  III. — Displaced  summation  for  the  four  quadrants,  II  and  III  being  turned  end  for  end,  and 
the  apex  of  the  curve  taken  for  origin  in  each  quadrant. 
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Some  interesting  conclusions  are  deducible  from  these  results.  One  is  the 
asymmetry  of  the  curve  representing  the  relation  between  the  proper  motions  and  the 
numbers.  From  the  adoption  of  the  apex  of  the  curve  as  the  zero  point,  it  follows 
that,  were  the  curves  symmetrical,  the  numbers  would  diminish  equally  in  each 
direction.  Inspection  of  the  table  shows  that,  up  to  4".o  on  each  side  of  the  assumed 
apex,  the  number  of  negative  proper  motions  is  the  greater.  This  shows  that  the 
adopted  apex,  which,  from  what  has  been  said,  was  only  roughly  assumed,  does  not 
correspond  to  the  true  one,  the  latter  being  somewhere  between  o".o  and  —  i".o.  But 
from  5".  outward  in  each  direction  there  is  a  constantly  increasing  excess  of  positive 
over  negative  proper  motions.  Of  motions  having  a  value  of  14". 6  and  upward,  171 
are  positive  and  only  58  negative. 

The  cause  of  this  asymmetry  is  so  obvious  as  scarcely  to  need  statement.  It 
arises  from  the  unequal  distances  of  the  stars,  and  the  consequent  inequality  in  the 
parallactic  motion  of  the  stars  of  any  one  class.  The  feature  in  question  may  enable  us 
to  form  a  rude  estimate  of  the  amount  of  this  inequality  of  distance  in  the  following 
way:  The  number  of  stars  whose  proper  motion  exceeds  H"-5  being  1 71,  we  find  in 
the  negative  part  of  the  table  a  proper  motion  —  /<0,  such  that  the  number  of  stars 
having  this  or  a  greater  proper  motion  shall  be  also  171.  We  find  by  count 
/x0  —  —  7". 3,  as  against  +  i4"-6  for  the  stars  in  the  opposite  direction.  Now,  it  may 
be  assumed  that  the  absolute  values  of  the  proper  motion  in  the  case  of  the  171  stars 
whose  motion  is  positive  is  the  same  as  in  the  case  of  the  1 7 1  stars  whose  motion  is 
negative.  But  the  difference  between  the  two  limiting  values  of  the  observed  proper 
motions  is  14". 6 —  7"-3  =  7"-3-  To  this  we  should  add  the  double  of  the  parallactic 
motion  for  the  apex  of  the  curve,  which  is  about  3".o.  This  gives  io".3  for  twice  the 
inferior  limit  of  the  parallactic  motion  of  the  342  stars  having  the  greatest  absolute 
proper  motion,  and  thus  5".  2  for  the  limit  of  the  parallactic  motion  itself.  Our  con- 
clusion is  that  out  of  the  2,874  stars  included  in  our  count,  about  342,  or  nearly  one- 
eighth,  have  parallactic  motions  of  5".  2  and  upward. 

But  it  must  be  noted  that  the  parallactic  motion  thus  found  is  not  the  total  motion, 
but  only  the  average  through  the  whole  circle  of  right  ascension  and  in  the  direction  of  the 
right  ascension.  The  maximum  parallactic  motion  in  right  ascension  for  the  same  system 
of  stars  is  found  by  multiplying  this  motion  by  n  2  =  1.57.  This  will  carry  the 
parallactic  motion  of  the  stars  whose  direction  makes  a  right  angle  with  the  direction 
of  the  apex  of  solar  motion  up  to  about  8".o.  It  should  be  still  further  multiplied  by 
sec.  D  to  get  the  total  absolute  parallactic  motion,  which  will  then  come  out  about  io".o. 
This  will  be  simply  the  inferior  limit  of  the  motion  for  the  group  of  342  stars  in  ques- 
tion, and  not  the  mean  for  all  of  them.  At  the  other  limit  the  motions  correspond, 
approximately,  to  the  apex  of  the  curve.  This  gives  us  for  the  more  distant  stars  of 
the  class  which  we  have  been  considering  parallactic  motions  of  i".5,  2".3,  and  3".o, 
respectively.  These  should  be  taken  as  the  motions  corresponding  to  the  great  mass 
of  the  more  distant  stars,  and  not  to  the  limit.  The  latter  is  not  capable  of  exact  dis- 
covery, because  a  few  scattering  stars  might  be  placed  at  a  great  distance  without 
materially  changing  the  result. 

It  is  not  claimed  that  the  preceding  conclusions  approximate  to  quantitative  pre- 
cision.   The  numbers  assigned  are  only  to  be  considered  as  slight  improvements  on 
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merely  qualitative  statements.  The  fact  that  the  apex  of  the  curve  does  not  approxi- 
mate to  a  zero  value  of  the  proper  motion,  but  to  an  effect  of  parallactic  motion 
amounting  to  nearly  3".o,  seems  to  show  with  some  conclusiveness  that  the  stars 
included  in  our  count  do  not  extend  out  indefinitely,  but  that  the  great  mass  of  them 
are  included  within  a  comparatively  limited  sphere — a  sphere  at  the  boundary  of  which 
the  parallactic  motion  would  seem  to  be  little,  if  any,  less  than  2".o. 

In  order  to  decide  upon  the  best  course  to  take  in  assigning  weights  as  a  function 
of  the  amount  of  the  proper  motion,  we  must  compare  the  law  of  distribution  with 
the  theoretical  distribution  according  to  the  normal  law  of  error.  A  very  slight  exam- 
ination will  show  that  if  we  endeavor  to  represent  the  actual  distribution  by  a  curve 
following  the  normal  law  the  closest  representation  will  leave  an  excess  of  actual  num- 
bers for  very  small  and  very  large  values  of  the  proper  motion,  while  the  intermediate 
numbers  will  be  deficient.  This  tendency  of  distribution  in  residual  deviations  is  so 
well  known  that  no  great  interest  will  attach  to  such  a  comparison.  What  I  have 
done  is  to  represent  the  curve  of  distribution  as  closely  as  possible  near  its  apex  by  a 
normal  curve,  and  then  note  the  deviations  of  the  two  curves  as  we  recede  from  the 
apex  in  each  direction.    The  result  is  shown  in  Table  IV. 

The  numbers  given  in  the  first  two  columns  of  this  table  are  taken  unchanged 
from  Table  III.  They  therefore  give,  for  each  value  of  a  proper  motion  in  the  direc- 
tion of  right  ascension  found  in  the  first  column,  the  number  of  stars  having  that 
value  of  the  proper  motion,  the  mean  effect  of  the  parallactic  motion  being  first  sub- 
tracted from  the  whole  series  in  the  way  already  described,  and  the  positive  direction 
being  always  taken  in  the  direction  of  the  parallactic  motion. 

An  examination  of  the  three  central  numbers  of  the  table  shows  that,  if  we 
represent  these  numbers  as  abscissa  of  the  curve  of  normal  error,  the  elements  of  this 
curve  will  be  based  on  a  probable  error  of  about  ±  1".  1 1,  and  the  apex  of  the  curve 
will  correspond  to  a  proper  motion  of  about  —  o".i5.  The  intervals  of  one  second, 
therefore,  correspond  to  0.90  of  the  probable  error.  Constructing  a  scale  of  abscissae 
in  which  the  probable  error  shall  be  regarded  as  unity,  the  base  of  the  first  ordinate 
on  the  negative  side  of  the  origin  being  placed  at  —  0.33,  and  the  first  one  on  the  posi- 
tive side  at +  0.5  7,  we  may  erect  successive  equidistant  ordinates  at  the  points  shown 
in  the  column  E  of  the  table. 

Next,  from  the  table  of  the  probability  of  errors,  such  IX  A  of  Chauvenet's 
Spherical  and  Practical  Astronomy,  Volume  II,  we  have  twice  the  probabilities  that  an 
error  will  be  contained  between  the  origin  and  each  of  the  abscissae  laid  down  in  the 
third  column.  These  probabilities  are  found  in  the  fourth  column,  P.  The  differences 
found  in  the  column  4  are,  therefore,  twice  the  probabilities  that  an  error  will  be  con- 
tained between  any  two  consecutive  points  on  the  axis. 

Multiplying  these  numbers  by  the  factor  694,  taken  arbitrarily,  so  as  to  satisfy 
the  condition  of  representing  the  three  central  numbers,  we  have  the  numbers  in  the 
next  column,  which  show  the  distribution  of  the  proper  motions  between  the  several 
limits  as  it  ought  to  be  according  to  the  normal  law  of  error  when  the  constants  are  so 
taken  as  to  nearly  satisfy  the  conditions  laid  down. 

Lastly,  the  deviations  from  the  actual  numbers  are  given. 
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Table  IV. — Comparison  of  actual  and  normal  distribution  of  the  excesses  of  the  proper  motions  of 
Bradley  stars,  in  the  direction  of  B.  A.,  over  the  mean  parallactic  motion  for  magnitudes  4.0  and 
upward. 
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The  value  of  the  probable  error  wliich  represents  that  portion  of  the  curve  near 
the  apex  is  thus  found  to  be 

Probable  error  =  ±  1".  1 1 . 

This  result  is  so  remarkable  a  one  that  I  should  hesitate  to  base  any  conclusions  upon 
it  were  it  not,  as  we  shall  soon  see,  completely  confirmed  by  the  distribution  of  the 
proper  motions  in  declination.  The  remarkable  feature  is  that  this  probable  error  is 
scarcely  greater  than  the  probable  deviations  of  the  numbers,  arising  from  the  variety 
of  values  of  the  parallactic  motion  which  have  been  combined,  and  from  the  probable 
errors  of  the  determination.  That  is  to  say,  had  1,388  stars  in  the  Bradley  catalogue 
no  proper  motion  except  the  parallactic  motion,  and  were  these  motions  determined 
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by  observation  and  statistically  tabulated  in  the  manner  we  have  done,  we  should 
expect  just  about  such  a  distribution  as  is  shown  in  this  normal  curve,  the  probable 
error  ±  i".ii  being  the  result  of  all  the  causes  just  mentioned.  If  we  should  now 
add  as  many  more  stars  having  all  values  of  actual  proper  motion  between  yet  wider 
limits  and,  to  fix  the  ideas,  suppose  them  equally  distributed  between  these  limits,  say 
between  +  1 5".  and  —  1 5".  per  century,  the  effect  would  be  to  rather  diminish  the 
sharpness  of  the  curve  at  the  apex,  and  therefore  to  increase  the  value  of  the  probable 
error  requisite  to  represent  the  central  portion  of  the  curve. 

We  thus  see  an  indication  that  about  one-  half  the  stars  in  Auwers-Bradley  of 
the  fourth  and  higher  magnitudes  have  no  appreciable  proper  motion,  except  what  is 
due  to  the  parallactic  motion.  If  they  have  any  absolute  proper  motion,  it  can  be 
only  a  fraction  of  a  second  per  century. 

The  only  grounds  I  see  on  which  this  conclusion  can  be  invalidated  is  that  there 
may  be  an  accidental  accumulation  of  proper  motions  near  the  value  o.  It  may  be 
remarked  that  in  Auwers-Bbadley  a  value  o  is  sometimes  assigned  to  a  proper  motion 
simply  because  the  data  were  insufficient  to  fix  upon  a  definite  positive  or  negative 
motion.  But  it  must  be  noted  that  all  these  cases  have  been  excluded  in  the  count. 
Moreover  the  values  which  lie  near  the  apex  of  the  curve  are  not  the  value  zero  of  the 
observed  proper  motions,  but  the  values  which  correspond  to  the  parallactic  motion.  It 
is  quite  true  that  in  a  chance  distribution  of  about  300  quantities  accidental  deviations 
of  a  magnitude  sufficient  to  change  the  form  of  the  curve  are  quite  possible.  The  con- 
clusion cited  must  not  therefore  be  regarded  as  a  proven  fact,  but  only  as  a  conclusion 
indicated  with  a  greater  or  less  degree  of  probability  by  statistics. 

Section  VIII. 

Statistical  distribution  of  the  proper  motions  in  declination. 

I  have  not  studied  the  proper  motions  in  declination  by  the  statistical  method  so 
fully  as  those  in  right  ascension,  because  it  is  not  so  easy  to  deduce  results  from  them. 
I  have,  however,  tabulated  the  distribution  of  the  proper  motions  in  declination  as 
given  by  Auwers  without  the  application  of  my  systematic  corrections.  As  the  uncor- 
rected values  will  answer  all  the  purposes  now  in  view,  the  results  and  comparison  are 
taken  from  a  paper  in  the  Astronomical  Journal  and  tabulated  below. 

As  to  magnitude,  the  stars  were  divided  into  four  groups  as  follows: 

Group  I,    Mag.  3.4  and  brighter ; 

"    n,   "  3.5  to  4.9; 

"      III,    "     5.0  to  5.9; 

"      IV,    "     6.0  and  upward. 

It  will  be  seen  that  there  is  an  excess  of  proper  motions  between  o".  and  —  1". 
per  century  defining  the  apex  of  the  curve,  and  as  well  marked  as  in  the  case  of  the 
right  ascensions.  The  negative  position  of  the  apex  is  evidently  due  to  the  parallactic 
motion  in  declination. 
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T^fiLB  V. — Distribution  of  Auwees's  centennial  proper  motions  in  declination,  according  to 
magnitude,  with  probable  distributions  according  to  the  normal  laic  of  error. 
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Section  IX. 

Formation  and  solution  of  equations  front  individual  stars. 

The  very  simple  formulae  for  the  equations  of  condition  between  the  four  unknown 
quantities  have  already  been  quoted.  Preliminary  to  the  formation  of  the  equations, 
all  the  stars  of  Auwers-Bradley  were  arranged  according  to  their  magnitudes  as  given 
in  Auwers's  work.  All  stars  up  to  magnitude  2.9  were  included  in  one  group,  because 
they  are  too  few  in  number  to  admit  of  any  discussion  of  smaller  groups.  Stars  of 
magnitude  3.0  to  3.9  were  tabulated  as  of  magnitude  3,  and  so  with  the  following  mag- 
nitudes up  to  7,  which  included  all  stars  of  the  7th  and  higher  magnitudes.  Separate 
normal  equations  were  then  formed  for  each  of  the  magnitudes. 

The  study  of  the  statistical  distribution  of  the  proper  motions  just  made  shows 
that  we  here  have  a  case  where  the  assignment  of  weights  as  a  function  of  the  deviation 
on  the  principles  already  developed  by  me  might  be  applicable.  I  have,  however,  felt 
that  not  enough  would  be  gained  by  the  application  to  justify  the  investigation  it  would 
involve.    I  have  therefore  adopted  the  usual  method  of  simply  excluding  all  proper 
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excluded  iu  investigations  on  the  subject.  In  fact,  all  the  stars  remaining  in  the  equa- 
tion will  be  among  those  hitherto  classified  as  having  small  proper  motions. 

The  limits  have  been  set  as  a  function  of  the  right  ascension  in  order  that  the 
effect  of  the  parallactic  motion  might  be  included  in  determining  them  in  each  sepa- 
rate case.  Properly  speaking,  the  limits  should  not  be  those  of  the  absolute  terms  in 
the  equation  of  condition,  or  the  absolute  proper  motion,  but  should  be  determined  by 
the  residuals  left  after  the  solution  of  the  equations  of  condition.  The  limits  thus 
derived  for  the  right  ascension  are  given  in  the  following  Table  VI: 


Table  VI. — Limits  outside  of  which  proper  motions  in  B.  A.  are  excluded. 


R.  A. 

Mag. 

1-2 

3 

4 

5-7 



h 

// 

// 

//  // 





//  // 

—  12  to  +24 

-  9+18 

—  7+15 

—  64-12 

\ 

— 12 

+  24 

—  9  18 

—  7+15 

—  6+12 

2 

—  13 

+  23 

— 10  17 

-  8+14 

—  7+11 

3 

-I4 

+23 

—11  16 

—  9+13 

-  7fii 

4 

-15 

-f-22 

—12  15 

—  9+12 

•--  8+11 

5 

—  l6 

+20 

—12  15 

— 10+ 1 1 

—  8+10 

6 

—  17 

+  19 

—13  H 

— 11-fii 

—  94-9 

7 

—  20 

fi6 

—14  «3 

— 12+10 

—104  8 

8 

—22 

+  15 

—  15  12 

—  124  9 

— 10+-  8 

9 

—22 

+  14 

— 17  11 

-13+  8 

—ii4-  7 

IO 

-23 

+  13 

—17  10 

-14+  8 

-124-  7 

ir 

-  24 

+  12 

—18  9 

-i44  8 

-124-  6 

12 

-24 

4  12 

—18  9 

-14+  7 

—12+  6 

13 

24 

+12 

—18  9 

—144  8 

—124-  6 

14 

-23 

+  13 

-17  10 

-134  8 

—11+  7 

15 

—22 

+  14 

—16  11 

-134-  9 

—"4-  7 

:6 

-21 

-M5 

—16  11 

-'34  9 

—Il4-  8 

17 

—20 

+  .6 

-15  12 

—12+  9 

-10+  8 

18 

—17 

+  '9 

-14  13 

— 12+ 10 

-94-9 

19 

—l6 

+21 

->3  14 

— 1 1  -f  10 

—  8+10 

20 

—5 

4-22 

—12  15 

—  IO-j-  II 

-  8+10 

21 

—14 

+23 

-11  16 

—  94  12 

—  7+11 

22 

—  '3 

+23 

— 10  17 

—  94  13 

—  74-12 

23 

— 12 

+24 

-  9  18 

-  8+14 

—  64-12 

In  the  case  of  the  declinations  the  limits  were  of  course  the  same,  except  that 
the  allowance  for  parallactic  motion  was  different  and  was  generally  smaller  than  in 
the  case  of  the  right  ascensions. 

The  general  rule  followed  was  to  exclude  a  star  if  its  proper  motion  in  either 
coordinate  was  outside  the  limits  of  the  table.  It  is,  indeed,  clear  enough  that  since 
the  proper  motion  has  no  relation  to  coordinates  the  exclusion  should  be  determined 
by  the  absolute  amount  of  the  proper  motion,  irrespective  of  the  question  whether  it 
chance  to  coincide  with  the  direction  of  one  of  the  coordinates.  For  this  reason  the 
rules  of  exclusion  were  not  entirely  absolute,  but  were  in  general  about  these :  Stars 
markedly  outside  the  limit  in  the  case  of  either  coordinate  were  always  excluded ;  if, 
however,  nearly  on  the  limit,  they  would  be  included  when  the  motion  in  the  other 
coordinate  was  small.  Moreover,  if  there  were  several  stars  in  any  region  of  the 
heavens  near  the  same  positive  or  negative  limit  of  exclusion,  some  would  be  included 
and  others  excluded.  When  the  motions  in  both  coordinates  approximated  to  the  limits, 
though  included  within  them,  the  star  was  excluded. 
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It  appears  by  count  that  the  number  of  excluded  stars  and  the  total  number  in 
the  several  classes  are  as  follows: 


Number  of  stars. 
Total.  Excluded. 


Mag.  1  and  2 
3 
4 
5 
6 

7 


87 
189 

43o 
949 
1.249 
277 

3.  18' 


23 
54 
'°3 
218 

215 
41 

654 


It  will  be  seen  that  a  little  more  than  one-fifth  the  entire  number  of  stars  are 
excluded. 

Section  X. 
Formation  of  the  normal  equations. 

In  forming  the  normal  equations  I  employed  in  most  cases  the  method  of  equiv- 
alent factors  proposed  by  Gr.  P.  Bond  in  the  Memoirs  of  the  American  Academy  of 
Arts  and  Sciences  about  the  middle  of  the  present  century,  but  rarely  applied.  It  is 
in  some  sort  an  intermediate  method  between  the  rigorous  treatment  by  least  squares 
and  the  system  sometimes  used  of  forming  normal  equations  from  conditional  equa- 
tions by  using  no  factors  except  +  i,  o,  or  —  i  to  form  each  normal  equation. 
Though  we  may  call  the  method  intermediate,  it  is  not  to  be  regarded  as  midway 
between  the  two  systems,  as  its  deviation  from  rigor  is  practically  inappreciable. 

The  method  amounts,  in  fact,  to  allowing  a  very  slight  change  in  the  weights  with 
which  the  several  conditional  equations  enter  into  the  normal  equations.  In  forming 
the  latter  each  conditional  equation  is  multiplied  by  the  product  of  the  coefficient  of 
the  corresponding  unknown  quantity  into  the  weight  of  the  equation.  Instead  of 
using  the  rigorous  value  of  this  factor  we  use  some  convenient  quantity  near  it,  say 
the  nearest  integer  or  the  nearest  tenth.  The  normal  equation  is  the  sum  of  the  prod- 
ucts of  the  conditional  equations  by  these  slightly  modified  factors.  In  this  method 
the  factors  [oft]  and  [6a]  are  not  rigorously  equal,  and  must  therefore  be  formed 
separately. 

In  the  present  case  I  have  taken  as  the  factors  ten  times  the  coefficient  of  each 
unknown  quantity  in  each  normal  equation.  Thus  we  have  the  following  normal 
equations  and  solutions  from  the  stars  of  each  class: 

Normal  equations  and  solutions  for  the  Bradley  stars. 


R.  A. 


402  ty  +  31  X+  l9Yrr+  753 

+  30    +350    —  77    =  +  156 
+  18    —  77    +301.4  =  —  1604 
Sp    =  +  2.2 
X    =— 1.0 


Dec. 


Mag.  1-2.9 


47«/>  +  8Y+  49X  +  nZ  =  +  10 
+  8     +80+21  o  542 

+  49  +21  +107  —  18  --—  118 
+  11  0—18  -L453  ^+1601 
Sp  =  —  0.8 
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.  Normal  equations  awl  solution*  for  the  Bradley  stars — Continued. 


R..A. 

// 

904^+iooX—  15Y--  +  278 

+ 100   +  850  —  46  =  —  98 

—  15    —46+501  —  —  2609 

Sp  =  +  0.27 

X  =  -0.43 

Y  =-5.21 


2165^/+     12X—  404  Y=+  1727 
+  18     +1822    +    39    ^-  378 
—  405     +    45    +1485    =  —  4658 
Sp    =  +  0.23 
X    =  —  o.  14 
Y  =-3.07 


4793^  +  472X-  653 Y  =  +  2880 


+  468  +3926 
—  660     +  270 


+  268 
+  3440 
Sp 
X 
Y 


7>55ty+  897X-  403 Y 
+  895     +5258    +  283 
—  410     +  260  +5144 
Sp 
X 
Y 


=  -  388 

=  -  9525 
=  +  0.22 
=  +  0.06 
=—2.73 


+  261 1 
^+  1555 
=  —  13552 
=  +  0.16 
=  +  0.42 
=  —  2.77 


16373/+  496X—  U4Y=^+  961 
+  492     +1310    +116  =+286 
—  119     +  »7    +1065    —  —  2451 
Sp    —  +  0.34 
X    =  +  0.29 
Y    =  —  2.30 


Mag.  3.0-3.9 


Dec. 


+  3 
+  60 
—  6 


78^+    3Y+  60X 
+  13 
+  138 
+  3i 


+  173 
+  13 
—  8 


+  31 

+  ,°# 
X 
Y 
Z 


6Z  =  —  93 
8   =-  743 
=  —  191 

=  +  2595 
=  +  0.50 
- 1.78 
=  —  404 
=  +  2.52 


Mag.  4-0-4-9 


24oty  +  10Y  +  129X  -  176Z 
+  12     +399-2    +  42 
+  127     —    1    +354    -  77 
—  171     +44—80  +25^5 

X 
Y 
Z 


~=+  39 
—  —  886 

=  -  94 
=  +  4291 
=  +  1-72 
=  —0.49 
=  —  2-47 
=  +1.84 


Mag-  S-o-5-9 


555«/>-  45Y+  472X-  318Z, 

—  54  +912  —  129  —  152  = 
+  458     —125    +1031     —  287  - 

—  308     —142    —  282    +5322  = 

dp  , 

X  - 

Y  = 

Z  = 


-+  334 

—  2448 

—  81 
+  8804 
+  2.00 
-0.83 

—  2.68 
=  + 1-58 


Mag.  6.0-6.9 


835s/—    3°Y+  667X—  211Z  — +  jioo 

—  36  +1236  —  92  —  678  =-  —  2798 
+  646     —    93    +1397    —  316  =+  752 

—  199     —  668    —  297    +7763  -=  +  12454 

6p  -  +  1.56 

X  —  +0.14 

Y  =  —  1 .01 

Z  =  +  1.56 


Mag.  7.0- 


tj2Sp—  11Y+160X—    58Z=  + 

—  11  +356  —  32 
+  154     —  24  +314 

—  59     —214    —  50 


,  200 

—    213      =r—  956 


—  53 
+  1709 

sp 

X 
Y 
Z 


=  +300 
=-  +  3°47 
—  +  1.06 

=  +  0.59 
=  —  1.62 

=  +  1.63 


Section  XI. 

Combination  of  the  partial  normal  equations  into  a  single  system. 

I  do  not  regard  the  results  of  the  jireceding  solutions  as  those  which  are  to  be 
finally  combined  in  a  definitive  result.  Separate  solutions  for  the  different  magnitudes 
have  been  made  in  order  to  ascertain  whether  there  is  any  systematic  difference 
between  the  precessional  motion  as  derived  from  stars  of  different  magnitudes,  or  as 
derived  from  right  ascensions  and  declinations,  respectively. 

There  is  evidently  a  marked  excess  in  the  value  of  Sp  derived  from  the  declina- 
tions over  that  derived  from  the  right  ascensions.    This  discrepancy  will  be  discussed 
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subsequently.  There  is  no  evidence  of  a  progressive  change  with  the  magnitude.  It 
is  true  that  the  results  from  the  few  brighter  stars  are  markedly  discordant;  but  this 
is  to  be  expected,  not  only  from  the  small  number  of  the  stars,  but  from  their  consid- 
erable miscellaneous  proper  motions.  When  the  results  from  the  right  ascensions  and 
from  the  declinations  are  combined,  it  will  be  seen  that  the  discordance  is  materially 
diminished. 

We  might  suspect  from  the  first  four  values  of  8p,  derived  from  the  declinations, 
that  there  is  a  progressive  increase  with  the  magnitude.  But  we  find  that  this  does 
not  continue  through  the  two  next  higher  magnitudes.  On  the  whole  we  can  not  say 
that  there  is  any  evidence  of  well-marked  progressive  change  with  magnitude.  We 
shall  therefore  in  the  subsequent  solutions  combine  all  the  magnitudes.  The  factors 
with  which  this  is  to  be  done  are  next  to  be  considered. 

What  we  should  really  do  is  to  multiply  the  conditional  equations  by  such  factors 
(square  roots  of  the  weights)  as  will  reduce  the  probable  errors  of  the  several 
equations  to  the  same  value.  But  these  probable  errors  are  nothing  more  than  the 
miscellaneous  deviations  of  the  proper  motions.  It  follows  that  if  we  form  a  series 
of  numbers  inversely  proportional  to  these  miscellaneous  residual  proper  motions  for 
the  several  magnitudes  of  the  stars,  such  numbers  will  form  the  appropriate  factors  of 
our  equations  when  the  latter  are  written  in  the  form 

«    .  sin  a  v    cos  a  v  » 
a  op  A  A  I  —  u  cos  o 

P  P 

,  _  .  cos  a.  sin  8  v  sin  a  sin  8  v  cos  6  „  , 
b  op  A  A  A  Y  Z—  u 

P  P  P 

We  have  still  to  consider  the  divisor  p  which  enters  into  the  equations  and 
expresses  the  distances  of  the  several  classes  of  stars.  There  being  no  method  of  deter- 
mining p,  we  are  obliged  to  assume  that  for  each  class  of  stars  it  is  inverselv  propor- 
tional to  the  mean  deviations  of  the  proper  motions,  and  therefore  directly  proportional 
to  the  factors  just  defined.    It  follows  that  if  we  write  the  equations  in  the  form 

a  p  8p  -f-  sin  a  X  —  cos  a  Y—  /x  p  cos  8  =  p  n 

b  p  Sp-\-  cos  a  sin  8  X  +  sin  a  sin  8  Y  —  cos  8  Z  —  fi'  p  —  p  n' 

the  equations  will  all  be  of  equal  weight.  We  must  therefore  determine  the  factors  p 
so  that  they  shall  be  approximately  the  inverse  of  the  proper  motions.  We  might 
equally  well  take  the  parallactic  motion  as  determining  p  for  each  class  of  stare. 
What  I  have  actually  done  is  to  find  the  deviations  in  various  ways,  and  to  consider 
also  the  conclusions  of  others  on  the  subject.  Fedorenko  (Astr.  Nach.  Vol.  XLV,  pp. 
81-84)  l°ug  since  found  that  the  proper  motions  of  stars  from  magnitude  1  to  8  were 
nearly  the  inverse  of  the  magnitudes.  This  approximate  law  was  confirmed  by 
Auwers  in  his  discussion  of  the  proper  motions  of  the  stars  contained  in  his  Berlin 
zone.  We  can  scarcely  regard  this  law  as  a  real  law  of  nature.  Its  most  obvious 
defect  is  that  it  would  make  the  proper  motions  of  stars  of  the  magnitude  o,  infinite. 
We  must  therefore  expect  that,  in  the  case  of  the  higher  magnitudes,  as  the  first  and 
second,  the  law  will  give  the  proper  motions  too  large.    Again,  in  the  case  of  the  very 
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faint  stars  the  proper  motions  included  are  not  the  mean  proper  motions  of  all  stars, 
but  only  of  a  selected  number.  When  all  stars  are  included  we  must  concede  that  the 
product  will  be  less  than  that  given  by  the  rule  as  fainter  stars  are  taken.  But  the  rule 
is  probably  very  nearly  true  for  the  magnitudes  with  which  we  are  concerned,  espe- 
cially those  from  the  fourth  to  the  seventh.  Comparing  the  rule  with  the  general 
probabilities  of  the  case,  with  the  average  deviations  shown  by  the  statistical  distribu- 
tion of  the  proper  motions  and  the  several  values  of  the  parallactic  motion,  I  have 
reached  the  conclusion  that  the  following  values  of  p  form  as  good  an  approximation 
to  the  required  factors  as  we  can  readily  obtain: 

Mags.  1-2.9;     9  =  0.4 

"  3-0-3-9;  =°-6 
"     4.0-4.9;  =0.8 

"  5-0-5-9;  =  1.0 
"  6.0-6.9;  =  1.2 
"     7.0-         =  1.4 

From  the  equations  last  written  it  will  be  seen  that  iu  forming  the  normal  equa- 
tions we  should  in  each  partial  normal  equation  already  given  multiply 

[a  a],  [a  n],  [b  b],  [b  n]    .    .    .    .    by  p2 

and  the  other  products  by  p,  or  by  unity,  according  to  the  degree  to  which  p  enters 
into  the  product  of  each  pair  of  coefficients  in  the  equations. 

When  the  coefficients  of  the  sets  of  normal  equations  found  on  pages  31-32  are 
multiplied  by  these  factors  and  combined  into  a  single  set,  we  have  the  following 
results.  The  right  ascensions  and  declinations  are  first  taken  separately;  then  they 
are  combined  by  addition: 

Combined  equations  from  all  Bradley  stars. 

R.  A.  " 

20079^)+  2324X—   i62iY=  +  9849 
23i7    +i35l6   +    583    =+  1579 

-  1645  +   570  +11937  =—35151 


Declination.  " 

2284<Sp—      82Y+1654X—      792Z  =  +'2302 

—  IOO      +3157    —    221     —    IOO9    =  —  8516 

+  1605    —  209  +334'    —    720  =+  1004 

—  764    —  988  —  696   +1882 1  =4-33644 

R.  A.  and  Dec.  combined. 

22363^4-  3978X—  1703Y—  792Z=+i2i5i 
3922    +16857    +    374   —    720  =+  2583 

—  '745    +    349    +15094   —  1009  =—43667 

—  764    —    696   —    988   + 1 882 1  =+33644 
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The  following  are  the  solutions  which  we  obtain  by  treating  the  three  sets  inde- 
pendently: 

R.  A.  Dec.  Combined. 

//  //  // 

Sp  =  +  0.25l  +I748  +0-357 

X=+o.203  —0.297  +0.203 

Y— — 2.920  — 2.108  — 2-745 

Z=  .     .  .         +1-737  +1.666 

As  there  can  be  but  one  value  of  the  several  quantities  X,  Y,  and  Z  for  any  given 
system  of  stars,  the  final  values  of  these  quantities  are  those  which  we  should  use  in 
obtaining  the  value  of  dp  from  each  separate  partial  normal  equatiou  in  Sp.  From  the 
factors  we  have  employed  it  will  be  seen  that  the  preceding  values  correspond  to  mag- 
nitude 5.5.  We  reduce  them  to  the  other  magnitudes  by  dividing  by  the  same  factors 
which  we  used  in  the  multiplication,  and  thus  obtain  the  following  results: 

Values  of  X,  Y,  and  Z  for  the  several  groups. 

Mag.    1-2.9  3.5  4.5  5.5  6.5  7.+ 

//  //  //  //  //  // 

X  +0.51  +0.34  +0.25  +0.203  +0.17  +0.14 
Y  —6.86  —458  —3-43  —2-745  —2.29  —1.96 
Z         +4.16      +278      +2.08      +1.666      +1.39  +1.19 

We  now  substitute  these  values  of  X,  Y,  Z  in  each  partial  normal  equation  in 
Sp,  whether  in  right  ascension,  declination,  or  the  combination  of  both,  in  the  combined 
equations  from  all  the  right  ascensions,  and  in  the  combined  equations  from  all  the 
declinations.  We  thus  have  the  -following  separate  results  for  dp  from  the  stars  of 
each  separate  magnitude  in  right  ascension  and  in  declination,  and  from  the  combina- 
tion of  both  motions,  together  with  the  final  results  from  all  the  right  ascensions,  all  the 
declinations,  and  the  combination  of  both  coordinates. 

Separate  and  combined  results  when  the  definitive  values  of  X,  Y,  and  Z  are  stibstituted 
in  the  normal  equations  in  Sp. 


R.  A. 

Dec. 

Both. 

n 

// 

// 

I  -  2 

Sp  =  +  2.16 

Sp  —  —  0.10 

<5p  —  +  1.92 

3 

+0.19 

-  1.05 

+  0.09 

4 

+  0.16 

+ 170 

+  0.19 

5 

+  0.21 

+ 1. 16 

+  0.31 

6 

+  0.22 

+  I-35 

+  0-35 

7 

+0.41 

+  1.30 

+  0-53 

All 

+0-245 

+  I-340 

+  0-357 

The  results  are  as  accordant  as  we  could  expect  when  derived  from  the  stars  of 
the  separate  magnitudes.  But  there  is  clearly  a  systematic  difference  between  the 
results  derived  from  the  right  ascensions  and  those  from  the  declinations.  In  the 
case  of  the  latter  the  result  is  the  more  noteworthy  from  its  discordance  with  that  of 
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L.  Stbuve,  who,  from  substantially  the  same  data,  obtained  a  result  which,  when  madf 
comparable  with  the  above,  would  imply  a  negative  correction  to  the  results  deriv  ed 
from  all  the  declinations.  As  a  prelude  to  considering  the  cause  of  the  discrepancy, 
and  to  deriving  the  most  probable  combination,  we  have  to  consider  the  possible  effect 
of  a  community  of  proper  motions  among  many  stars  in  different  regions  of  the 
heavens.  In  connection  with  this  the  consistency  between  the  results  derived  from 
the  stars  in  different  directions  will  have  to  be  considered. 

Section  XII. 

Classification  of  the  stars  by  regions  and  zones. 

The  results  obtained  in  the  preceding  section  have  been  derived  by  giving  each 
individual  star  equal  weight.  But  it  is  well-known  that  a  community  of  proper 
motions  frequently  exists  among  small  groups  of  stars — the  Pleiades,  for  example.  It 
is  evident  that  in  such  cases  the  entire  group  should  receive  only  the  weight  of  a 
single  star,  or  perhaps  two.  Were  we  to  seek  out  the  special  cases  of  this  kind,  so 
far  as  they  could  be  determined  from  the  equations,  they  would  be  so  few  as  not  to 
materiallv  affect  our  result.  What  I  have  therefore  done  is  to  divide  the  heavens  into 
regions,  each  forming  a  trapezium  of  one  or  two  hours  of  right  ascension  and  extend- 
ing through  15°  or  30°  of  declination,  and  find  the  mean  proper  motion  of  all  the 
stars  in  each  trapezium.  A  process  somewhat  similar  to  this  was  followed  by  L. 
Struve  in  combining  his  equations  from  the  separate  stare.  From  the  mean  of  all 
the  stars  in  each  trapezium  he  formed  a  single  equation  of  condition,  assuming  the 
mean  proper  motion  to  correspond  to  the  mean  position.  He  then  gave  to  each  equa- 
tion a  result  proportional  (making  abstraction  of  differences  of  magnitude)  to  the 
number  of  stai-s  which  the  trapezium  contained.  Thus,  his  results  are  the  same  as  if 
he  had  found  and  solved  the  conditional  equations  for  each  star  separately,  as  I  have 
done  in  the  preceding  section.  I  consider  the  right  ascensions  and  declinations 
separately. 

Section  XIII. 

discussion  of  the  right  ascensions. 

When  equations  are  combined,  on  the  supposition  of  a  possible  community  of 
proper  motions  within  each  region,  the  question  of  combining  stars  of  different  mag- 
nitudes offers  some  difficulty.  Such  a  common  motion  presupposes  that  the  stars  are 
nearly  at  the  same  distance,  and  therefore  have  approximately  the  same  proper  motion. 
It  would  seem,  therefore,  that  in  such  a  case  all  stars  should  be  given  equal  weight, 
irrespective  of  their  magnitude.  But  it  is  certain  that  there  are  many  scattered  stars 
to  which  the  hypothesis  would  not  apply.  What  I  have,  therefore,  done  in  the  case  of 
the  R.  A.'s  is  to  confine  the  discussion  to  stars  of  magnitudes  4.0  and  upward.  Equal 
weight  was  then  given  to  each  star,  so  that  the  mean  proper  motion  within  each  zone 
and  region  is  that  of  the  stars  in  question  within  the  region.  The  zones  are  150  wide, 
except  in  the  case  of  that  south  of  —  1 50.  Here  I  have  included  all  stars  as  if  belong- 
ing to  a  single  zone.  The  mean  proper  motions  thus  derived  from  the  second  mem- 
bers of  the  equations  of  condition  are  written  in  the  table  below.    For  example,  in 
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Table  VII. — Equations  for  precession  from  proper  motion  in  R.  A.  by  regions,  etc. — Continued. 


Zone+150  to  +30° 

Zone  +-300  to  +450 

R.  A. 

6 

=  +22^° 

1  =  37^° 

o 



P 

X 

y 

n 

Wt. 

A 

X 

y 

// 

Wt. 

77* 

.86 

+•13 

—■99 

-- 

+3-°«3 

•7 

76 

-I--I3 

-■99 

— 

+  0.5.3 

6 

11  \C 
**7t 

+  .38 

-.92 

— 

2.8* 

/ 

.82 

+  .38 

-.92 

I.4.. 

6 

3/ 72, 

■94 

+  .61 

—•79 

- 

2.3* 

8 

.88 

+.6i 

-.79 

— = 

2.9„ 

6 

Z.1  \£ 

527i 

•97 

+  ■79 

-  .61 

■■- 

1-835 

8 

.92 

+•79 

—.61 

= 

i. a, 

c 

°77* 

■99 

+■92 

-.38 

■- 

4-i35 

X 

■V5 

+  .92 

-■38 

= 

1-4.5 

Hi  V. 

-t-99 

—.13 

- 

i-9*. 

8 

0*7 
•y/ 

+  -99 

-•13 

o.8l6 

6 

at  \L 

1  00 

+  ■99 

+  13 

— 

+o-5j6 

Q7 

•y/ 

+  ■99 

o-5„ 

6 

1  *  72 

+  .92 

+.38 

- 

8 

•yo 

+  .92 

+.38 

— 

+  o.9,„ 

c 

l27% 

•97 

+■•79 

+.6r 

-  2.4j, 

8 

•y-* 

+•79 

+.61 

—  °-9m 

6 

IA1  vz 

•94 

+.61 

-K79 

— 2.4.4 

.88 

+  .61 

+•79 

-  o.8u 

6 

TC*T  \C 

l577* 

•9° 

+.38 

+.92 

—3-2.6 

.82 

+  .38 

+.92 

-  3-3= 

7 

I  "71  \£ 

.86 

+-«3 

+  •99 

—1. 9-» 

6 

76 

+-13 

+  99 

—  5-°» 

2 

18*7  \L 

.82 

-13 

+  ■99 

—2.9,, 

8 

70 

—•'3 

+■99 

—  2.7„ 

c 
0 

202% 

•79 

-•38 

+  •92 

-2.5.0 

5 

•63 

—38 

+  92 

-  7-84 

3 

217  \k 

76 

-  .61 

-•79 

—I. I8 

C 

0 

-0" 

-.61 

+■79 

—  10.8, 

1 

2327a 

•  72 

—■79 

+  .61 

-  3-6.8 

7 

-54 

-■79 

+.61 

-  1.8,, 

6 

247M 

•71 

-.92 

+  •38 

-1-4.5 

6 

•50 

■92 

f.38 

—  1-9.5 

6 

26272 

.70 

-  -99 

+  13 

+0.2,5 

6 

■48 

—99 

+■13 

-  i-3« 

5 

277  yi 

.70 

■99 

-•13 

o.8„ 

6 

.48 

—•99 

-13 

+  o.4,9 

7 

292  % 

71 

-.92 

-.38 

«-3» 

8 

•50 

-.92 

-.38 

+  °-9m 

6 

3°7# 

•72 

-•79 

-  .61 

0.6,, 

7 

•54 

~l9 

-.61 

2-Qi, 

7 

322  % 

•  76 

-.61 

-•79 

3-4., 

6 

•58 

—.61 

-•79 

I.O.. 

6 

337H 

•79 

-•38 

—.92 

0.7, 

4 

-63 

-•38 

-.92 

3-4.3 

6 

352X 

.82 

-•13 

-•99 

+2.5,6 

6 



.70 

—  13 

—•99 

+  3-5.. 

6 

Zone  +450  to  +60° 

Zone  +600  to  +900 

R.  A. 



l-52#° 

c  —  700 



0 

P 

X 

y 

// 

—  - 
Wt. 



A 
P 

X 

y 

// 

Wt. 

in 

.60 

+■'3 

—■99 

-H.6,3 

g 

•3° 

+■13 

—•99 

+  2.3„ 

7 

'68 

+  .38 

—.92 

2.2,5 

6 

aft 

+.38 

-.92 

2.6,5 

7 

•/4 

+  61 

-■79 

1.6,0 

6 

Cyl 
■34 

-r.6i 

-•79 

0.1,8 

7 

.80 

+■79 

—.61 

3-°.. 

g 

6l 

+•79 

-.61 

2-34 

3 

°7/2 

8c 
■°5 

+■92 

-•38 

I.O„ 

.66 

+  92 

-.38 

+O.I, 

2 

■°7 

+•99 

-13 

+2.0,, 

7 

.68 

+-99 

—  13 

—O.49 

5 

0*7  \L 

87 

.07 

+  •99 

+•13 

-0.9,0 

r 
0 

.68 

■\  -99 

+.13 

 0.6,o 

c 

0 

1 w/2 

8c 

•°5 

+  ■92 

+  .38 

-  0.1,, 

5 

.66 

+.92 

+.38 

+O.I„ 

5 

127/2 

.80 

+  •79 

+.61 

+3-5s 

4 

.61 

+•79 

+.61 

-0.5, 

5 

lA272 

•  74 

+.61 

+•79 

+0.36 

4 

•54 

+  .61 

+  -79 

—4-28 

5 

TC7  l«£ 

•57/2 

.68 

+•38 

+■92 

-1.66 

4 

.40 

+.38 

+.92 

— 2.2. 

4 

+  •13 

+•99 

-•-3? 

4 

■3° 

+•13 

+•99 

—0.2, 

187^ 
10/72 

4 .99 

— O.I, 

c 
0 

27 

-•'3 

+•99 

—  I-9u 

20272 

•44 

-:3 

f  .92 

-2. 18 

5 

•7 

--38 

+.92 

O.O, 

1 

217/2 

.^6 

—.61 

4  -79 

2-95 

4 

.09 

-.61 

+•79 

-2.8, 

2 

23272 

•3' 

•79 

i  .61 

f6.8, 

1 

-f.02 

-•79 

4.61 

+2.4. 

2 

247K 

.26 

-.92 

+.38 

-2.4* 

4 

-•°3 

-.92 

+.38 

— I.lj 

4 

"  262  'A 

•24 

-•99 

+  ••3 

-  1.05 

4 

-.06 

-■99 

+  -«3 

-2.1, 

2 

■24 

-•99 

-••3 

+0.6, 

5 

-  .06 

—■99 

— 0.2„ 

6 

292^ 

.26 

-.92 

-•38 

+O.Is 

5 

-.03 

.92 

-.11 

-+0.48 

5 

307  M 

•31 

-•79 

—.61 

1.6,5 

6 

—  .02 

-•79 

-.61 

1.1,0 

7 

32272 

•36 

-.61 

-•79 

2.7,, 

6 

.09 

—.61 

—79 

1.2,8 

7 

337  % 

•44 

-•38 

-.92 

1.6,, 

6 

•17 

-■38 

-■92 

1-9.5 

6 

352# 

•52 

—■13 

—•99 

-t-2.4.7 

7 

+  .27 

—•13 

-•99 

+1.6* 

7 

By  the  third  process  we  should  find  a  result  substantially  the  same  as  that  already 
derived  by  the  solutions  from  the  individual  stars.  The  deviations  which  would  arise 
from  the  center  of  the  trapezium  not  always  coinciding  accurately  with  the  center  of 
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gravity  of  the  stars  and  from  the  different  magnitudes  having  equal  weight  would  be 
likely  to  produce  only  a  small  change  in  the  result 

In  order  to  bring  out  the  differences  sought  for,  I  have  made  a  solution  for  each 
zone  separately. 

In  the  solution  based  on  equal  weights  for  each  trapezium  I  have  omitted  the 
extreme  southern  zone,  owing  to  the  irregular  distribution  of  the  stars,  and  also  the 
northern  one,  owing  to  its  small  weight.    To  eliminate  the  parallactic  motion  I  have 

simply  added  the  equations  in  each  of  the  other  zones  with  the  results: 

// 

Zone  — 15°  to       o°;    2i.SSp  =  —  0.2 

o°to  +  i5°;    21.8         +  7.1 

+  '5°  to  +30°;    20.3         +  34 

+  30°  to +45°;  174  +  3-3 
+  45°to+6o°;    13.4  +  18.1 

From  these  we  derive  the  normal  equation 

1850  S  p  =  +520". 

Whence, 

dp  —  +o".28. 

The  second  or  adjusted  system  of  weights  gives  rise  to  the  following  normal 
equations: 

// 

598.7  <*i>+  58.1X  —  33.7Y  =  +  225.0 
476.6    +  17.9   =+  76.1 
456.8   =  —  1231.2 

These  give 

// 

dp  =  +0.20 
X  =  +O.24 
Y  —  —  2.69 

This  I  regard  as  the  preferable  result  of  the  method,  and  have  therefore  made  a 
comparison  of  the  partial  results  from  the  seven  separate  zones  of  Dec.  This  is  done 
by  taking  the  partial  normal  equation  in  p  given  by  each  zone,  and  substituting  the 
concluded  values  of  X  and  Y  so  as  to  derive  a  value  of  6 p.    The  results  are  as  follows: 


Zones. 

Equations. 

>i 

-37° 

to  -25° 

93  4  SP 

36.4  +  30.8X4  i4Y  = 

+ 

40.O;  Sp 

=  + 

0.43 

-15° 

to  o° 

123.2 

26.4+  2.8 

+  5-9  = 

+ 

I  1.2 

+ 

0.09 

o° 

to  4  15° 

1 2 1 .4 

51-9-  8.3 

+  4-9  = 

+ 

36.7 

+ 

0.30 

+  15° 

to  +  300 

1 21.0 

38.0  —  21.5 

+  1.9  = 

+ 

27-7 

+ 

0.23 

+  30° 

to  +  45° 

74.1 

U0.3  —  1 7.7 

+  8.3  = 

l6.2 

0.22 

+  45° 

to  +  6o° 

46.3 

60.2  —  23.2 

+  7-9  = 

+ 

33-3 

+ 

0.72 

+  6o° 

to  +  8o° 

19-3 

1.8  —  21.0 

4  34  = 

12.3 

0.64 

Total, 

598.7  SP 

225.0  —  58.1 

+  337  = 

+ 120.4;  SP 

=  +  0.201 
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By  the  third  system,  in  which  the  weights  are  proportional  to  the  number  of  stars, 
we  have  the  following  normal  equations  and  solutions: 

// 

1568^+  187X—   i68Y=+  796 
1209    —   17  =+231 
1086    =  —  2,975 

// 

X  =  +  O.I2 
Y  =  -2.7I 

p  =  +  0.204 

This  is  less  by  o".04  than  the  result  from  each  star  separately,  a  difference  which 
arises  partly  from  the  fact  that  in  this  solution  stars  of  all  magnitudes  have  equal 
weights.  If  we  combine  the  four  results  on  page  3  5,  for  the  magnitudes  4  to  7,  namely, 
".16,  ".21,  ".22,  ".41,  giving  each  a  weight  proportional  to  the  number  of  stars,  we 
have  o".225.    We  thus  have  the  following  three  sets  of  values  for  p  by  the  three 


separate  systems  of  treatment: 

// 

All  equal  weights,  jp  =  +  o.28 

Adjusted  weights,   +0.20 

Separate  stars,   +  0.24 


These  differences  are  so  slight  as  to  show  that  the  result  will  not  be  materially 
altered  by  allowing  for  a  supposed  community  of  proper  motions  in  different  regions. 
The  slight  deviations  from  the  former  result  seem  rather  to  increase  the  deviation  from 
the  result  of  the  declinations. 

A  singular  feature  is,  however,  shown  by  the  results  when  we  arrange  them  by 
right  ascensions.  It  is  not  possible  to  give  a  definite  result  for  the  precession  from 
each  separate  hour  of  right  ascension,  owing  to  the  effect  of  the  parallactic  motion. 
But  from  the  mean  result  of  any  two  opposite  hours  of  right  ascension  the  parallactic 
motion  may  be  eliminated,  always  supposing  that  the  stars  of  equal  apparent  magni- 
tudes are  at  equal  distances.  I  have  proceeded  in  the  following  way  to  make  a  com- 
parison of  mean  proper  motions  derived  from  regions  of  equal  and  opposite  values  of 
the  parallactic  motion.  Four  regions  of  right  ascension  are  included  as  follows: 
Taking  the  right  ascension  A  of  the  solar  apex  as  1 8  h.  30  m.  the  adopted  regions  are 
lunes  designated  as  (a)  (b)  (c)  and  (d). 

(a)  Right  ascensions  within  2  J  hours  ±  from  A  +  900; 

(b)  "  "  "      "      "  A+1800; 

(c)  "  "  "      "      "  A  —  900; 

(d)  "  "  "     "     "  A. 

It  will  be  seen  that  each  of  these  regions  includes  five  hours  of  right  ascension, 
leaving  a  space  of  one  hour  between  adjacent  regions.  Within  the  four  regions  the 
mean  proper  motions  are  taken,  assigning  to  those  within  each  trapezium  the  weights 
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employed  in  the  preceding  equations  of  condition.  The  results  for  the  mean  proper 
motions  in  the  direction  of  R.  A.  are: 


Zone. 

(a.) 

(b.) 

(c) 

(d.) 

#(a+c) 

#(b+d) 

Diff. 

ft 

// 

// 

tr 

// 

// 

->n° 
—  20 

t  3-35 

,1      r  nQ 

—  3-49 

1    r\  88 

  O.O/ 

-r  °-95 

+  1-°5 

-15° 

o° 

+  2.38 

—  O.3O 

-3-4° 

+  142 

-O.SI 

+  0.56 

+  1.07 

o°  + 

15° 

+  2.21 

+  0.43 

—  2.27 

+  1. 16 

—  0.03 

+  0.79 

+  0.82 

30° 

+  2.40 

+  O.60 

-2-39 

+  0.37 

—  O.OO 

+  0.48 

+  0.48 

45° 

+  2.34 

+  0  53 

-4.49 

+  O.I4 

—  I.O8 

+  0.34 

+  142 

45° 

6o° 

+  i-9° 

+  O.96 

-  i-55 

—  0.02 

+  O.I8 

+  0.47 

+  °-29 

Mean 

+  2.43 

+  o-55 

—  2.93 

+  O.66 

—  0.25 

+  0.60 

+  0.85 

I  have  stopped  at  +  60  0  of  declination  because  no  result  of  value  for  our  purpose 
could  be  derived  beyond  that  limit.  Were  the  stars  affected  only  by  a  parallactic 
motion  equal  in  the  general  mean  in  the  case  of  different  magnitudes  the  two  columns 
before  the  last  should  (leaving  out  the  effect  of  differences  in  the  coefficients  of  the 
correction  to  the  precession)  be  equal  in  the  general  mean.  But  there  is  evidently  a 
difference  of  a  systematic  character,  the  mean  proper  motions  near  the  meridians  which 
include  the  solar  apex  being  larger  than  those  forming  a  right  angle  with  that  region. 
One  way  of  accounting  for  this  would  be  by  supposing  that,  taking  stars  of  equal 
apparent  magnitude,  those  in  region  (a)  which  lies  on  the  two  sides  of  the  vernal 
equinox  are,  in  the  general  mean,  farther  away  than  those  lying  on  the  two  sides  of 
the  meridian  of  twelve  hours  of  right  ascension.  This  amounts  to  the  same  thing  as 
supposing  that  stars  in  the  vernal  region  are  somewhat  larger  in  absolute  magnitude 
than  those  in  the  autumnal  region.  The  result  would  be  that,  comparing  equal  mag- 
nitudes, the  stars  (c)  will  have  larger  parallactic  motions  than  the  stars  (a).  The 
mere  probable  effect  of  the  accidental  inequalities  of  distribution  and  magnitude  would 
result  in  a  certain  inequality  of  this  kind,  but  the  uniformity  of  the  discrepancy  pre- 
vents our  attributing  the  difference  to  this  cause  alone. 

We  may  more  plausibly  attribute  the  difference  to  possible  systematic  errors  of 
Bradley's  right  ascensions,  those  near  o*  and  1 2*  of  R.  A.  being  too  large  relative  to 
those  near  6*  and  18*.  The  amount  of  the  correction  which  would  annul  the  discrep- 
ancy is,  approximately, 

4a  —  —  o'.<>33  cos  ia. 

That  so  large  an  error  as  this  could  have  affected  Bradley's  observations  certainly 
seems  unlikely.  Yet,  in  view  of  its  approach  to  constancy  through  8o°  of  declination, 
as  shown  by  the  column  of  differences,  it  seems  difficult  to  attribute  the  discrepancy 
to  any  other  cause.  That  systematic  inequalities  in  the  proper  motions  themselves 
could  attain  this  magnitude  and  follow  this  law  seems  very  unlikely. 

It  will,  however,  be  of  interest  to  examine  more  closely  in  what  way  the  discrep- 
ancy is  related  to  the  right  ascensions.  I  have  therefore  divided  the  sphere  into  eight 
lunes  of  right  ascension,  so  situated  that  the  central  meridian  on  the  initial  lune  should 
coincide  with  the  assumed  meridian  of  the  solar  apex,  in  eighteen  hours  and  thirty 
minutes  of  right  ascension.    The  normal  equations  for  p  formed  from  the  numbers 
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in  Table  VII  were  added  for  each  of  these  lunes  through  all  seven  zones,  with  the  fol- 
lowing results: 

(a)  (b) 


R.  A. 

//  // 

18.5; 

65.6  S  j) 

—  76.6X 

-11.6Y 

=  +  77-3 

2>  =  +  0.96;  -f-  1.1 1 

21-5; 

7i-5 

-  53-2 

—  68.8 

=  +197-3 

+  0.47;  +0.58 

0.5; 

744 

+  12.3 

—  90.4 

—  4-  224.2 

—  0.08;  —  O.IO 

3-5; 

88.7 

+  78.4 

-60.8 

=  +  201.5 

4-0.38;  +0.25 

6.5; 

100.5 

+  108.9 

+  13-4 

=  +  34-i 

+  0.47;  +0.32 

9-5; 

78.9 

+  53-7 

+  7I-3 

=  -163.3 

+  0.07:  —0.03 

12.5; 

59-8 

-  9-5 

+  71-7 

=  -214.7 

—  0.54;  —0.51 

15-5; 

59-3 

-  55-9 

+  4I-5 

=  —  132.0 

—  0.28;  —0.15 

598.7 

+  58.1 

-33-7 

224.4 

Except  on  the  meridian  and  autimeridian  of  the  solar  apex  the  several  values  of 
6p  to  be  derived  from  these  individual  lunes  will  depend  on  the  adopted  value  of  the 
parallactic  motions  X  and  Y.  If  we  use  the  values  of  X  and  Y  as  already  derived 
from  the  conditional  equations  we  should  take  the  values  for  mag.  6.0,  which  is  nearly 
the  mean  magnitude  of  the  stars  we  have  used.    We  shall  then  have 

// 

X  =  4-o.i9 

Y  =  -2.52 

and,  for  S  j>,  the  results  given  above  as  (a).  But  if  we  take  the  solar  apex  as  situate 
in  1 8  h.  30  m.  of  right  ascension,  we  shall  have 

X  =  Y  cot.  2  7 70. 5 

and  should  therefore  adopt 

// 

X  =  +  o.33 

Y  =  -2.52 

We  shall  then  have  the  numbers  given  in  column  (1)). 

We  may  also  combine  the  results  in  opposite  hours  of  right  ascension,  so  as  to 
very  nearly  eliminate  the  effect  of  the  parallactic  motion  from  each  pair  of  opposite 
lunes.  This  is  done  by  multiplying  either  equation  of  a  pair  by  such  a  factor  that 
when  added  to  the  opposite  equations,  the  coefficients  of  X  and  Y  shall  be  so  small 
that  the  uncertainty  of  the  parallactic  motion  will  not  materially  affect  the  result. 
Practically,  however,  it  will  suffice  to  take  the  mean  result  of  each  pair  of  opposite 
lunes.    We  thus  have  the  following  four  values  of  Sji: 

R.  A. 


h. 

h. 

// 

6.5  and 

18.5 

dp  —  -f-  0.72 

9-5 

21-5 

+  0.27 

12-5 

0.5 

-0.31 

15-5 

3-5 

+  0.05 

The  dependence  of  the  result  upon  the  region  of  right  ascension  from  which  the 
elements  of  the  problem  are  derived  is  brought  out  yet  more  strongly  by  this  com- 
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parison.  We  have  four  pairs  of  lunes — one  we  may  call  the  apical  pair;  that  at  right 
angles  to  it,  the  equinoctial  pair;  the  other  two,  intermediate  pairs.    We  have 

// 

from  the  apical  pair  8 p  —  -\-o.j2 

•"      "   intermediate  pairs         +0.1 6 
"      "   equinoctial  pair  —0.31 

If  we  take  a  mean  of  these  three  results,  using  the  weights  as  derived  from  the  normal 
equations,  the  weights  will  be  nearly  in  the  ratio  11,  20,  and  9,  and  the  result 

6p  =  4-o".2i, 

substantially  that  already  derived.  But  it  is  clear  that  since  the  result  of  the  apical 
pair  is  free  from  all  the  uncertainty  arising  from  the  parallactic  motion,  it  should 
receive  a  larger  proportional  weight.  How  much  larger  can  not  be  determined  by  any 
mathematical  process;  I  would  suggest  the  factors  1.0,  0.7,  and  0.5  as  those  which 
commend  themselves  to  my  judgment  The  weights  will  then  become  1 1,14  and  4.5, 
or  22,  28  and  9.    Tho  result  for^>  will  then  be 

dp  -  +  o".30 

If  we  judge  the  probable  error  of  this  result  by  the  discordances  just  shown,  it 
would  be  at  least  ±  o".2  5.  But  if  we  consider  the  motion  of  the  equinox  as  denned 
by  the  mean  of  the  Bradley  stars  to  be  correct  and  the  discordance  to  arise  from  sys- 
tematic errors  in  the  R.  A.'s,  then  the  discordance  will  cease  to  be  an  index  to  a 
probable  error,  and  the  simple  mean  +  o".2i  will  be  the  required  result.  But,  as  I 
have  said,  the  inequality  of  parallactic  motion  may  possibly  be  a  factor  in  the  result. 
I  therefore  consider  tho  result  o".30  to  be  that  given  by  the  li.  A.'s  of  the  Bradley 
stars. 

The  question  may  arise  whether  the  decisively  larger  value  given  by  the  apical 
pair  of  lunes  is  related  to  the  parallactic  motion,  or  grows  out  of  the  fact  that  a  large 
proportion  of  the  stars  in  these  lunes  is  situated  in  tho  Milky  Way.  This  ques- 
tion demands  a  special  study  of  the  galactic  stars,  with  a  view  of  determining  whether 
they  give  results  systematically  different  from  those  derived  from  other  stars. 

Section  XIV. 

Special  examination  of  galactic  stars. 

It  is  well  known  that  the  condensation  of  stars  in  the  region  of  the  Milky  Way 
begins  with  the  lower  magnitudes,  and  continually  increases  as  we  ascend  to  higher 
and  higher  magnitudes.  The  researches  of  Kapteyn  indicate  that  tho  galactic  stars 
are,  in  the  general  average,  farther  than  the  others ;  in  fact,  that  the  galaxy  may  be 
considered  as  an  annular  collection  situate  outside  the  general  lx>undary  of  the  rest  of 
the  system.  But  it  does  not  follow  from  this  view  that  any  large  proportion  of  the 
galactic  stars  in  the  Bradley  catalogue  really  belong  to  the  galaxy  as  thus  defined. 
Referring  to  Kapteyn's  researches  for  a  more  exact  statement  of  his  conclusions,  we 
shall  here  see  what  light  may  be  thrown  on  the  subject  by  an  examination  of  Auwers's 
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Berlin  A.  Gr.  catalogue  for  the  zone  +  150  to  +  200.  This  catalogue  is  especially 
valuable  for  our  purpose,  owing  to  the  very  exhaustive  discussion  of  the  proper 
motions  within  its  limits,  which  Auwers  has  there  given.  Within  this  zone  I  take  as 
the  galactic  regions, 

h.  h. 

Region  A,  5.0  to  8.0; 
Region  B,  18.0  to  21.0. 

These  two  regions  include  six  hours  of  right  ascension,  or  one-fourth  of  the  entire 
surface  of  the  zone.  The  corresponding  numbers  of  stars  contained  in  the  catalogue 
are, 

Stars  in  region  A,  1,797 
Stars  in  region  B,  i,a84 

Total  number  galactic  stars,  3,781 
Outside  of  the  galaxy,  6,008 

Total  number  of  stars,  9,789 

It  thus  appears  that  0.386  of  the  stars  of  the  catalogue  are  contained  with  the 
galactic  regions,  the  density  being  somewhat  greater  in  region  B. 

The  corresponding  statistics  for  the  number  of  stars  which  Auwers  has  found  to 
be  affected  by  an  appreciable  proper  motion  are, 

Proper  motions  in  region  A,  155 
Proper  motions  in  region  B,  202 

Total  in  the  galactic  regions,  357 
Outside  of  the  galaxy,  901 

Total,  1,258 

An  equal  distribution  of  the  1,258  proper  motions  would  give  157  within  each  of 
the  regions  A  and  B.  The  actual  number  in  region  A  falls  short  of  this  by  2,  and  that 
in  region  B  exceeds  it  by  45. 

On  the  other  hand,  were  the  numbers  proportional  to  the  stars,  we  should  have, 

In  region  A,  230  stars  with  proper  motion; 
In  region  B,  240  stars  with  proper  motion. 

It  will  be  seen  that  there  is  no  excess  in  the  number  of  proper  motions  in  region 
A.  In  the  case  of  region  B  there  is  an  excess  which  may  lead  us  to  suspect  some 
condensation,  yet  the  condensation  is  much  less  than  proportional  to  the  number  of 
stars.  The  conclusion  we  may  draw  is  that,  assuming  that  the  stars  having  sensible 
proper  motions  are  all  included  within  a  sphere  of  nearly  definite  radius,  the  stars 
within  this  sphere  do  not  show  any  well-marked  condensation  in  the  galactic  regions, 
unless,  perhaps,  in  region  B  near  1 9  h.  of  right  ascension,  and  therefore  do  not  belong 
to  the  galaxy  proper,  as  defined  by  Kapteyn. 
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the  value  of  q1  is  larger  for  the  higher  magnitudes  than  for  the  lower  ones.  Although, 
when  all  the  stars  are  considered,  this  parallactic  motion  is  somewhat  smaller  than  that 
already  derived,  the  difference  is  not  sufficient  to  show  any  marked  diminution  for  the 
case  of  galactic  stars. 

We  now  pass  to  the  more  southern  regions  of  the  galaxy  on  each  side  of  the 
sphere.  The  limits  of  these  regions  I  have  taken  from  Heiss's  Atlas  within  each  zone. 
By  adding  the  normal  equations  in  dp  for  each  region,  as  given  in  Table  VII,  we  have 
the  following  results: 


Zone. 

limits  of  R.  A 

to 

45° 

3 

h—  6h. 

,  15.2fy-1-14.5X 

-6.6Y 

=+17.2 

18 

—  22 

7.2 

-1 1-3 

-6.6 

=  +15-2 

-15° 

to 

30° 

5 

-  7 

16.0 

+  15-8 

+  0.0 

=+19.2 

18 

—  21 

10.6 

-134 

-5-7 

=  +13.8 

o° 

to 

15° 

5 

-  8 

18.6 

+  18.6 

+  1.9 

=  +  7-7 

17 

—  20 

13-4 

-14.9 

—  2.2 

=  +18.7 

-15° 

to 

o° 

6 

-  8 

8.2 

+  9i 

+  2.3 

=-  7-7 

17 

—  20  ' 

16.4 

-  16.6 

-30 

=  4-21.1 

-35° 

to  - 

■i5° 

7 

-  8  ; 

3-3 

+  40 

+  1.7 

=4-  0.4 

17 

-19  ; 

14.0 

-13-8 

+0.0 

=4-  5-6 

The  sum  of  these  equations  give  the  following  normal  equation  in  dp : 

1 22.9^=  11 1".2  +  8.0X4-  18.2Y 
or 

Sp  =  o".905  4-  0.065X  4-  o.  1 48Y 

Taking,  as  for  mag.  6, 

X  =  +  o.i7 
Y  —  —  2.29 

this  will  give 

Sp  =  +  o".58 

while  the  value  will  be  somewhat  larger  if  we  take  the  value  of  q  cos  D  found  from 
the  stare  of  the  galactic  region  alone. 

If,  to  eliminate  as  far  as  possible  the  parallactic  motion,  we  confine  ourselves  to 
the  galactic  stars  south  of  300  Dec,  we  shall  have 

100.5SP  -  78".8  +  1 1.2X  4  5.0Y  =  69".3 

and 

<5p  =  4-o".69 

On  the  whole  it  does  not  seem  that  we  are  here  dealing  with  anything  else  than 
the  systematic  changes  of  correction  already  found  for  the  meridian  and  antimeridian 
of  the  solar  apex.  This  periodic  difference  may  bo  connected  with  the  Milky  Way  or 
with  the  solar  motion,  it  is  impossible  to  say  which ;  all  that  we  can  do  at  present  is 
to  refer  it  to  positions  in  the  heavens. 
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Section  XV. 

Discussion  of  the  proper  motions  in  declinatian. 

The  results  of  the  preceding  discussion  of  the  right  ascensions  seem  to  indicate 
that  we  have  to  fear  systematic  errors  of  the  observations  and  data  more  than  irregu- 
larities of  proper  motion.  Hence,  in  testing  the  declinations  in  the  same  way  that  we 
have  the  right  ascensions,  I  first  took  zones  300  wide  instead  of  150,  and  in  each  zone 
belts  two  hours  wide  in  right  ascension.  I  also  varied  a  little  the  manner  of  taking  the 
means  within  each  trapezium ;  instead  of  an  indiscriminate  mean  of  all  the  stars  from 
the  fourth  to  the  seventh  magnitude,  I  assigned  to  the  mean  result  for  the  stars  of  each 
magnitude  a  weight  found  by  multiplying  the  number  of  stars  by  the  following  factors: 

Mag.,  4567 
Factor,  0.2     0.3     0.4  0.5 

The  results  of  the  process  are  shown  in  Table  VIII.  The  four  zones,  whose  limits 
are  given  at  the  top  of  the  table,  are  designated  by  the  letters  A,  B,  C,  and  D.  Here 
in  each  trapezium,  after  each  of  the  four  magnitudes,  is  given  first  the  mean  proper 
motion  in  declination  of  the  stars  of  that  magnitude,  followed  by  the  number  in  sub- 
script figures  of  stars  on  which  each  mean  depends.  This  is  followed  by  the  corre- 
sponding weight,  and  the  weighted  mean  is  taken  as  the  final  result  within  the 
trapezium. 

In  order  to  find  the  effect  of  the  systematic  corrections  which  I  have  applied  to 
Auwers's  proper  motions  in  declination,  the  same  means  were  taken,  using  Auwees's 
proper  motions  unchanged.    The  final  results  only  of  this  process  are  given. 

1  have  used  these  results,  not  to  obtain  a  separate  definitive  value  of  the  pre- 
cession, but  to  judge  of  the  probable  magnitude  of  the  systematic  errors  and  the 
effect  of  the  systematic  corrections  to  Auwers's  proper  motions.  Owing  to  the  com- 
binations of  all  the  proper  motions  within  each  region,  it  is  not  possible  to  make  any 
distinction  between  the  different  values  of  p  for  the  different  magnitudes.  Indeed 
from  the  way  in  which  the  stars  are  distributed  it  is  probable  that  the  effect  of  the 
different  values  of  p  will  be  but  slight.  Regarding  the  value  of  p  as  unity,  each 
proper  motion  gives  an  equation  of  condition  of  the  form 

cos  a.  (  An  +  X  sin  8)  -\-  Y  sin  a  sin  d  —  Z  cos  S  —  ft'. 

In  the  combined  equation  we  should  take  within  each  trapezium  the  mean  values 
of  the  sines  and  cosines  of  the  right  ascension  and  declination  which  enter  into  the 
expression.  I  have,  however,  considered  that  it  would  answer  our  present  purpose 
sufficiently  well  to  regard  these  mean  values  as  those  corresponding  approximately  to 
the  middle  point  of  each  trapezium.  In  fact,  the  general  agreement  among  the  num- 
bers seems  to  show  that  the  results  would  be  changed  very  slightly  by  any  probable 
deviation  of  these  values.  An  exception  may  possibly  be  made  in  the  case  of  the 
region  within  300  of  the  pole. 

Multiplying  each  of  the  twelve  equations  of  the  form  last  written  by  cos  a  and 
adding  the  products,  we  shall  have 

6  {A  n  +  X  sin  S)  —  2  pi'  cos  a. 
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The  results  for  the  four  zones  are  as  follows.  In  deriving  the  second  column  of 
numbers  from  the  first  I  have  taken 

X  =  +  0".20. 

It  must  be  noted  that  the  J'n  derived  from  Auwers's  equations  applies  to  the  Struve- 
Peters's  value  of  n. 


N'S 

AUWBRS'S 

Zone. 

A  n  +  X  sin  S 

An 

Wt. 

A'n  +  X  sin  D 

A'n 

tt 

11 

// 

11 

o° 

to  -36° 

+  0.59 

+  O.63 

2 

—  O.I3  — 

O.O9 

o° 

to +30° 

+  0.45 

+  O.4O 

4 

—  O.54  - 

o-59 

+300 

to  +6o° 

+  O.47 

+  0.34 

3 

—  O.26  — 

0.39 

+  6o° 

to  +90° 

+  O.81 

+  O.63 

2 

—  O.4O  — 

0.58 

Mean 

A  n 

=  +  0.47 

An  —  — 

■0.44 

or  8  p 

=  +  I.I7 

4  »  =  +  0.41 

S  p  —  +  1.02 

Applying  these  corrections  to  the  respective  values  of «,  we  find  the  following  mean 
results  for  the  value  of  this  quantity  at  the  epoch  1850 : 

// 

From  separate  stars,  100  n  —  2005.32 

"  zones  and  regions,  2005.26 
"     Auwers's  prop,  motions,  2005.20 

L.  Struve's  result  is  2004.55 

This  result  may  be  subject  to  some  doubt,  owing  to  the  breadth  of  the  zones  and 
the  variation  of  the  parallactic  motion  within  that  breadth.  Still  greater  doubt 
may  arise  from  the  fact  that  the  effect  of  the  parallactic  motion  could  not  be  rigor- 
ously taken  account  of  in  determining  whether  a  star  of  considerable  proper  motion 
should  or  should  not  be  included.  Since  the  limits  of  inclusion  are  properly  deter- 
mined only  by  the  magnitudes  of  the  residuals  after  all  corrections  for  precession  and 
parallactic  motion  are  made,  it  follows  that  a  star  which  should  be  excluded  when 
judged  by  the  absolute  amount  of  its  proper  motion,  might  be  included  when  the 
residual  proper  motion  became  known  and  vice  versa. 

A  yet  narrower  extension  of  the  limits  is  suggested  by  the  Btrength  of  the  evi- 
dence, already  presented,  that  nearly  one-half  the  Bradley  stars  have  no  motus 
peculiaris  in  one  coordinate  exceeding  i".o  per  century.  Could  we  with  certainty 
pick  out  these  stars  from  others,  it  would  be  well  to  depend  upon  them  exclusively  for 
the  determination  of  the  precession.  But  as  we  can  not  do  so,  an  extension  of  the 
limits  of  proper  motions  to  be  included  must  be  made,  and  we  can  then  find  no  stop- 
ping place  inside  db  l"  ■  In  order  to  use  proper  motions  as  near  as  possible  to  the 
final  residuals  in  determining  whether  a  star  should  or  should  not  be  included,  and 
in  order  to  test  the  whole  process,  I  proceeded  as  follows  with  a  new  and  nearly 
independent  determination  by  zones  150  wide  of  declination  and  single  hours  of 
right  ascension.  I  began  by  applying  to  Auwers's  proper  motions  in  declination  a 
A  P,  VOL  vih,  pt  1  4 
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precessional  correction  corresponding  to  a  diminution  of  Struve's  ioo  »  by  o".53,  and 
subtracting  a  parallactic  motion 

—  2". 50  sin  d  sin  a  4-  o".25  sin  S  cos  a  —  i"-50  cos  S. 

The  sphere  was  divided  into  zones  150  wide  north  of  the  equator,  while  south  of 
the  equator  only  the  single  zone  extending  from  o°  to  —  200  was  used.  Stars  brighter 
than  magnitude  4.0  were  omitted,  and  all  stars  above  this  magnitude  were  included 
with  equal  weight  except  so  far  as  exclusion  on  account  of  large  proper  motion  was 
concerned,  which,  no  distinction  of  magnitude  being  made,  led  to  the  exclusion  of  a 
larger  proportion  of  the  brighter  stars. 

As  to  the  effect  of  the  terms  of  parallactic  motion,  it  is  to  be  considered  that  the 
coefficient  of  cos  a  is  important,  because  it  is  in  each  zone  united  with  the  coefficient 
expressing  the  correction  to  the  precession.  Fortunately,  the  close  approximation  of 
the  solar  apex  to  eighteen  hours  of  right  ascension  makes  the  term  of  parallactic 
motion  depending  on  cos  a  quite  small;  we  may  consider  o".20  and  o".^o  as  its  prob- 
able limits.  Moreover,  it  is  diminished  by  being  multiplied  by  sin  S,  so  that  an  error 
in  its  assumed  magnitude  diminishes  to  zero  at  the  equator  and  changes  its  sign  south 
of  the  equator.  The  term  in  sin  a  is  without  effect,  except  so  far  as  the  stars  may  be 
unequally  distributed  near  the  two  sides  of  the  equinox.  The  adopted  coefficient  is 
so  near  that  obtained  for  the  mean  of  the  included  stars,  and  the  sum  total  of  the 
asymmetry  on  the  two  sides  of  the  equinox  is  so  small,  that  nothing  would  be  gained 
by  attempting  to  correct  this  motion  from  the  equations  in  question.  We  shall,  in 
fact,  be  nearer  the  truth  by  assuming  the  coefficient  to  be  a  known  quantity.  These 
considerations  have  led  me  to  apply  throughout  each  zone  of  declination  the  paral- 
lactic motion  corresponding  to  the  mean  declination  of  the  zone,  as  found  by  substi- 
tuting the  mean  value  of  sin  6  and  cos  8  for  each  zone  in  the  above  expression. 

Within  each  zone  the  parallactic  motion  in  declination  will  be  merged  with  the 
mean  systematic  correction  to  Auwers's  proper  motions  within  the  zone.  This  correc- 
tion is  derived  from  an  unpublished  paper  on  the  standard  system  of  declinations. 

The  following  are  the  expressions  for  the  corrections  actually  applied  to  the  proper 
motions  of  each  zone,  the  last  number  including  the  systematic  correction  to  the  A.  G. 
system  as  well  1".$  cos  S: 


0  0 

Zone  A;  —  20  to  0; 

+  0.57  cos 

a  —  0.40  sin  a  +  0.75 

.     B;  oto+15; 

+0.50  > 

+  0.20  ) 

+  0.32 

+  0.60 

•     C;  +  15  to +  30; 

+  0.44  ? 
+  0.04  i 

+  0.94 

+  050 

D;  +30  to +45;  +0.38  > 
—  0.42  s 

+  1.50 

+  0.20 

E;  +45  to  +  60; 

+  °-33 

+  1-97 

+  104 

F;  +6oto+75; 

+  0.30 

+  2.30 

+  0.93 

•    G;  +  75  to+  90; 

+  0.28 

+  2-45 

+  0.80 

I  began  by  applying  to  the  proper  motions  the  additional  terms  4S„,  derived 
from  a  comparison  of  the  Boss  standard  stars  witli  Auwebs-Bradley,  and  heretofore 
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considered  as  corrections  applicable  to  the  Bradley  declinations.  South  of  the  equator 
this  correction  is  zero.  North  of  the  equator  the  mean  value  within  each  zone  is  written 
under  the  coefficient  of  cos  a  up  to  45 °.  For  example,  in  zone  B,  the  correction  for 
precession  and  parallactic  motion  is  +o"-50  cos  a  and  that  for  ddm  is  -f-o".20  cos  a. 
North  of  450  45.  is  omitted.  After  including  it  in  the  first  three  zones  north  of  the 
equator  I  began  to  doubt  its  reality,  and  did  not  apply  it  in  the  remaining  zones. 
There  is  thus  a  slight  lack  of  similarity  in  the  results  actually  used,  but  as  the  effect 
of  the  change  in  each  zone  is  rigorously  to  change  the  result  in  the  processional  motion 
by  the  amount  of  the  coefficient,  no  confusion  need  result. 

The  proper  motions  having  thus  been  arranged  by  zones,  and  those  of  each  hour 
brought  together,  it  became  necessary  to  determine,  in  the  case  of  each  star,  whether 
it  should  or  should  not  be  included  on  account  of  large  proper  motion.  The  adopted 
rule  was  to  include  all  whose  residual  proper  motion  was  now  less  than  7",  and  to 
exclude  all  that  exceed  10".  Between  these  limits  a  star  was  retained  or  omitted 
according  to  whether  it  had  a  small  or  a  large  proper  motion  in  right  ascension,  and  also 
according  to  the  number  of  stars  having  motions  of  nearly  the  same  amount  When, 
as  frequently  happened,  two  or  three  neighboring  stars  within  the  same  zone  had  con- 
siderable proper  motions  of  nearly  the  same  magnitude,  some  would  be  struck  out  and 
others  retained.  In  the  case  of  the  Pleiades  the  common  proper  motion  of  the  group 
was  considered  as  having  the  weight  of  two  stars  only.  The  whole  process  was  con- 
ducted with  great  care  to  have  no  bias,  which  might  be  produced  by  any  knowledge  of 
the  effect  of  an  exclusion  upon  the  final  result. 

Two  other  points  in  which  the  process  differs  from  the  previous  one  were  that  only 
stars  within  four  hours  of  the  equinox  were  included,  it  being  considered  that  the 
weight  of  the  precessional  motion  derived  from  stars  within  300  of  a  solstice  was  too 
small  to  warrant  their  inclusion.  Yet  another  point  was  that  the  proper  motions  of  those 
stars  whose  declinations  were  not  determined  by  Bradley  were  left  out 

Table  IX. — Mean  residual  proper  motions  in  declination  for  Bradley  stars  of  magnitude  4.0  and 
upward,  after  subtracting  precession  (20".0511oos  a  for  1850)  and  a  parallactic  motion  of 
—2". 50  sin  S  sin  a+0".25  cos  a  —  l"J)  cos  S. 


Zone  A ;  — 200  to  o°. 

Zone  B;  o°  to  +150. 

ZoneC; 

+  i5°to+30°. 

ZoneD;  +300  to +450. 

R.  A. 

Mo 

No.  of 

Wt. 

Ho 

No.  of 

Wt. 

M« 

No.  of 

Wt. 

Mo 

No.  of 

Wt. 

h. 

// 

stars. 

n 

stars. 

// 

stars. 

// 

stars. 

20.5 

+0.63 

32 

11 

+O.I3 

22 

10 

—0.90 

19 

9 

—O.40 

22 

10 

21-5 

+  1. 12 

30 

11 

+O.40 

16 

8 

+0.1 1 

14 

8 

—  1. 15 

11 

7 

22.5 

+0.64 

36 

12 

+  I.07 

15 

8 

-0.15 

11 

7 

—0.68 

'3 

8 

23-5 

+0.40 

28 

11 

+O.I  I 

19 

9 

—0.97 

18 

9 

+0.04 

11 

7 

o.5 

-0.77 

17 

9 

+0.43 

29 

11 

+0.21 

23 

10 

—O.77 

>4 

8 

1.5 

— 0.06 

17 

9 

—0.67 

15 

8 

— 0.62 

27 

11 

—O.08 

11 

7 

2.5 

—1.28 

19 

9 

—0.41 

12 

7 

+0.15 

27 

11 

—115 

13 

8 

3-5 

+0.54 

19 

9 

+0.48 

11 

7 

-1.84 

20 

10 

— O.71 

9 

6 

8.5 

— 0.60 

16 

8 

-(-0.26 

16  • 

8 

+0.23 

33 

11 

+0.34 

14 

8 

9-5 

-+0.20 

25 

10 

+0.13 

21 

10 

—0.84 

17 

9 

+0.2I 

11 

7 

10.5 

4-I.20 

19 

9 

+0.23 

26 

11 

— 0.06 

20 

10 

+O.I2 

21 

10 

11.5 

+  I.60 

10 

7 

—0.58 

20 

10 

— 0.41 

»3 

8 

—  I.20 

3 

12.5 

—O.81 

15 

8 

—I-3I 

14 

8 

—0.95 

29 

11 

+O.OI 

9 

0 

13-5 

+O.30 

28 

11 

— 2.  II 

8 

6 

-1.63 

11 

7 

-(-O.02 

4 

3 

14.5 

-I.84 

21 

10 

—1-37 

7 

5 

—I.OI 

10 

7 

+3.90 

1 

1 

15-5 

—  I. OO 

26 

11 

—  I.OI 

17 

9 

— 1.51 

16 

8 

—O.86 

10 

7 
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Table  IX. — Mean  residual  proper  motions  in  declination  for  Bradley  stars  of  magnitude  4.0  and 
upward,  after  subtracting  precession  (20".0511cos  a  for  1850)  and  a  parallactic  motion  of 
sin  6  sin  a  +  0".25  cos  a  — 1".5  cos  8— Continued. 


 1 — 

Zone  E;  +45°  to  +  6o° 

ZoneF; 

+60°  to +750. 

ZoneG; 

+75010+90°. 

R.  A. 

Mo 

No.  of  Wt. 

M« 

No.  of 

Wt. 

f*o 

No.  of 

Wt. 

h. 

tt 

stars. 

rr 

stars. 

11 

stars. 

20.5 

— o.  19 

16 

8 

+0.24 

13 

7 

—  2.24 

4 

3 

21-5 

+0.27 

11 

7 

— 0.76 

7 

5 

-(-0.03 

5 

4 

22.5 

—0.40 

12 

7 

— 0. 14 

8 

6 

-f-0.29 

8 

6 

23-5 

+0.66 

18 

9 

— 0.43 

17 

9 

—0.40 

4 

3 

o-5 

+0.24 

12 

7 

— 0. 12 

12 

7 

— 0. 15 

10 

7 

i-5 

+0.72 

12 

7 

+O.68 

20 

9 

—0.24 

6 

5 

2-5 

+0.46 

13 

7 

+  «-32 

12 

7 

+  I-57 

4 

3 

3-5 

-1.28 

12 

7 

+O.62 

4 

3 

8-5 

+  1.34 

5 

4 

—O.81 

9 

6 

9-5 

— O.II 

8 

6 

+  I.OI 

8 

6 

—1.77 

2 

2 

10.5 

+0.14 

5 

4 

—2.97 

3 

3 

+o-45 

5 

4 

11.5 

+0.58 

4 

3 

+4.20 

1 

+1.80 

12.5 

—0.32 

9 

6 

-0.85 

6 

5 

5 

4 

13-5 

—0.87 

11 

7 

—1.80 

1 

+  1.58 

14-5 

—0.08 

5 

4 

+1.65 

2 

2 

2 

2 

15-5 

—0.20 

2 

2 

—0.89 

2 

2 

TABLE  X. — Coefficients  for  solutions  by  least  squares. 


Zone. 

[«0 

[ztfcos  a] 

[if  cos3  a] 

O  /*o] 

[W  UtiCOS  <t] 

// 

A 

155 

+  6.5 

+108.6 

+  3-9 

+16.8 

B 

135 

+  1.1 

96-3 

—25-1 

+45.8 

C 

146 

+  3-9 

105.2 

-93-8 

+18.4 

D 

106 

+  13-8 

73- 1 

-37-7 

-37-3 

E 

93 

+  196 

67.8 

+  3-2 

+  1-7 

F 

79 

+24.4 

56.6 

+  2.9 

+  12.3 

G 

45 

+  '5-3 

35-o 

+  3-4 

-  9-2 

The  mean  results  within  each  zone  and  each  hour  of  right  ascension  are  shown  in 
Table  IX.  It  is  then  assumed  that  the  mean  of  the  proper  motions  corresponds  to  the 
central  point  of  each  trapezium,  which  is  in  the  middle  of  each  hour  of  R  A. 

The  second  column  of  each  block  of  the  table  contains  the  mean  of  the  proper 
motions  after  the  exclusions  and  corrections  already  described.    This  is  followed  by 
the  number  of  included  stars  and  the  corresponding  weight,  the  latter  being  deter 
mined  on  the  hypothesis  that  the  stars  in  each  trapezium  may  have  a  proper  motion 
equal  to  one-fourth  the  mean  proper  motion  of  the  individual  stars. 

From  each  proper  motion  in  each  zone  an  equation  of  condition  is  formed  between 
the  correction  to  the  precession  in  declination  and  the  constant  correction  required  to 
all  the  declinations  within  the  zone.    This  equation  is : 

x  +  cos  a4n  —  fiQ 

The  reason  for  omitting  the  term  in  sin  a  has  already  been  given. 
From  the  equations  of  each  zone  normal  equations  in  x  and  J  n  are  formed.  The 
coefficients  of  these  equations  are  shown  in  Table  X.    The  equations  were  solved 
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assumed  that  the  mean  of  these  motions  should  vanish  when  account  is  taken  of  the 
parallactic  motion.  It  will  be  seen  that,  within  the  zone  i5°-300  there  is  strong 
evidence  of  a  systematic  error;  this,  however,  must  be  discussed  elsewhere. 

Section  XVI. 

Comparison  of  the  results  of  proper  motions  in  declination  with  those  of  L.  Struve. 

The  most  remarkable  feature  of  the  results  so  far  obtained  is  the  discordance 
between  the  values  of  the  precessional  motion  derived  from  the  right  ascensions  and 
from  the  declinations.  I  have  already  remarked  that  this  discrepancy  extends  through 
all  the  magnitudes,  from  the  fourth  upward.  It  is  not  materially  diminished  by  any 
mode  of  treatment  of  either  the  right  ascensions  or  the  declinations.  When  a  result 
from  right  ascensions  is  derived  by  using  the  Lalande  stars,  the  discrepancy  is  dimin- 
ished. But  I  do  not  consider  that  the  result  from  these  stars  has  more  than  a  small 
fraction  of  the  weight  to  be  attributed  to  the  Bradley  stars.  We  can  not  attribute  the 
difference  to  the  variation  of  personal  equation  with  magnitude,  because  we  have 
reason  to  suppose  that -this  error  is  one  which  affects  the  modem  more  than  the  older 
observations,  and  which  would,  therefore,  make  the  precessional  motion  derived  from 
right  ascensions  come  out  too  large.  What  adds  to  the  embarrassment  is  that  from 
substantially  the  same  material — namely,  a  comparison  of  Auwers-Bbadley  with 
modern  observations — L.  Stbuve  derived  a  value  of  n  so  much  smaller  than  mine  that, 
accepting  his  result,  the  discordance  would  be  reversed.  In  fact,  the  result  of  the 
direct  solution  of  the  equations  derived  from  declinations,  as  found  on  page  1 8  of  his 
paper  *  is, 

4n  —  —  i  ".09. 

My  provisional  motion  of  n  is  less  than  the  Strove-Peters  value  used  by  L.  Struve 
by  o".85.  Hence,  L.  Stbuve's  result,  considered  as  a  correction  to  my  provisional 
value,  would  be 

// 

z/n  —  —  0.24 
8  p  — —  0.60 

This  value  of  p  diverges  farther  in  the  one  direction  than  my  own  does  in  the  other. 
The  difference  between  our  values  of  n  is 

N.-L.  S.  =  +  o".68. 

This  enormous  discrepancy  seems  to  show  that  one  result  or  the  other  must  be  seriously 
in  error,  either  in  the  numbers  used  or  in  the  mode  of  treatment.  I  have,  therefore, 
sought  to  investigate  its  source,  and  applied  every  test  I  could  think  of  to  my  own 
methods  and  conclusions. 

The  first  question  that  would  arise  is  whether  the  discrepancy  extends  through 
all  the  zones  of  declination,  or  is  peculiar  to  some  of  them.  I  therefore  made  a  rough 
solution  of  Stbuve's  normal  equations  for  /in  in  his  various  zones,  using  his  concluded 

♦Menioires  de  L'Academie  Imperiale  des  Sciences  de  Saint-Petersbourg.  Septieme  serie.  Tome  XXXV. 
Saint-Petersbourg,  1887. 
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parallactic  motion,  and  assigning  equal  weights  to  each  of  his  trapezia.  The  assign- 
ment of  equal  weights  prevents  the  numbers  from  being  definitive  ones ;  but  the  errors 
thus  arising  are  of  no  importance  for  the  purpose  of  the  present  comparison,  which 
shows  that  the  discrepancy  extends  through  all  the  zones  of  declination,  unless,  pos- 
sibly, those  near  the  pole. 

The  stars  of  different  magnitudes  are  very  unequally  weighted  in  the  two  investi- 
gations, Stbuve's  diminishing  much  more  rapidly  with  the  brighter  magnitudes  than 
mine  do.  If  the  discrepancy  arose  from  this  cause,  my  results  from  stars  of  the  higher 
magnitudes,  say  6  and  7,  should  approach  more  nearly  to  those  of  Stbuve.  But  the 
actual  divergence  is  in  the  opposite  direction.  The  stars  of  magnitude  1  -  3  give 
negative  corrections  in  the  preceding  discussion.  Hence,  were  a  smaller  weight  given 
to  these  stars,  my  result  would  be  increased,  which  would  make  it  more  divergent 
from  the  other. 

I  have  applied  the  corrections  d  <Ja  to  Auwebs's  proper  motions.  These  make  the 
value  of  very  slightly  larger.  1  was  finally  led  to  reject  them  for  reasons  men- 
tioned in  the  preceding  section,  but  doing  so  does  not  remove  the  discordance. 

A  very  probable  cause  of  the  discrepancy  would  be  the  very  different  systems 
which  we  used  in  rejecting  stars  on  account  of  large  proper  motion.  Stbuve  rejected 
only  seven  stare  for  this  cause,  and  n  Cassiopeiae  is  not  included  in  the  list  of 
those  rejected.  It  might  therefore  well  happen  that  the  stars  which  I  have  rejected, 
but  which  he  included,  would  materially  affect  the  result.  With  the  view  of  deter- 
mining whether  this  was  the  cause  of  the  discordance,  I  have  formed  normal  equations 
in  the  four  unknown  quantities  from  those  stare  which  were  omitted  by  me,  but 
included  by  Stbuve.    The  results  and  solution  are  as  follows: 

Normal  equations  and  solutions  from  large  proper  motions  in  declination  used  by  L.  Struve, 

but  omitted  in  the  preceding  investigation. 

Mag.  4. 

QO.i<5p  —  2.3Y  +  60.6X  —  31.6Z  =—1856 
—  2.3  +125.0  —  3.5  +  32.0  —1253 
+  60.6  —  3.5  +127.3  —  24-6  -H45 
+  31.6     +  32.0    —  24.6    +7674  +8070 

Mag.  5. 

i88.q<Jp+  6.1Y  +  110.9X—  38.2Z=+  1181 
+  4.3  +261.5  -  77  -  50.5  -  3177 
+  106.5     —    7.2    +250.5    +    26.0       +  175 

—  34.2     -  48.8    +  27.9    +1662.0       +  17758 

Mag.  6. 

187.3$^—  10.7Y  +  121.6X  —     67.iZ=+  390 

—  10.8  +217.4  —  20.0  —  80.4  —  4966 
+  121.6     —  20.0    +283.7    —     98-7       +  l9°6 

—  67.0     —  80.4    —  98.6    +  1629.1  +11024 
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Normal  equations  and  solutions  from  large  motions  in  declination  used  by  L.  Struve,  but 
omitted  in  the  preceding  investigation — Continued. 

Mag.  7. 

38-3*P-  34Y+304X-  39-6Z=-  619 

—  34  +37-8  -14.8  -  564  —  321 
+  304     —148    4-557    —  21.8       —  282 

—  39-6     —564    —21.8    +297.7  +I951 

Combined  system. 

59i.3<Sp  —  13.3Y +347.9X  —  I994Z=—  658 
—  15-3  +6417  —  46.0  -  155.3  -10587 
+  343-5  -  45-5  +717  2  -  1191  +  '024 
-195-3     -!53-6    -1 171    +4356.2  +40174 

Solution. 

dp  =  +  0.52 
X  =  +  1.72 
Y  =  -  14.25 
Z  =  +  8.78 

There  is,  as  we  might  anticipate,  a  good  deal  of  discordance  between  the  results 
from  stars  of  different  magnitudes.  What  is  especially  curious  is  the  large  negative 
absolute  term  given  by  stars  of  magnitude  4.  This  arises  from  a  large  number  of 
stars  of  this  magnitude  around  the  vernal  equinox  having  a  negative  proper  motion. 
I  regard  this  accumulation  of  negative  proper  motions  as  purely  accidental. 

It  will  be  seen  that  the  value  of  p  derived  from  stars  of  large  proper  motion, 
while  less  than  that  from  stars  of  smaller  proper  motion,  is  yet  markedly  larger  than 
Stbuve's  result.  The  combination  of  these  stars  with  those  of  small  proper  motions, 
which  I  have  used,  would  diminish  the  value  of  n  by  only  a  few  hundredths  of  a 
second.    It  would  seem  that  the  discordance  can  not  be  accounted  for  in  this  way. 

It  would  seem  therefore  that  there  must  have  been  some  discordance  between  the 
proper  motions  in  declinations  found  by  L.  Struve  from  his  comparison  of  Auwebs- 
Bbadley  with  the  Pulkowa  declinations  for  1855,  and  those  derived  by  Auweks  from 
a  comparison  of  the  Bradley  catalogue  with  a  modern  system  for  1865.  The  method 
of  reduction  and  comparison  adopted  by  the  two  investigators,  as  described  by  them, 
was  substantially  identical.  What  we  have  therefore  to  do  is  to  find  whether  there 
can  be  any  systematic  error  in  Auwers's  proper  motions  as  found  in  his  Bradley  cata- 
logue. The  only  error  of  the  kind  which  would  materially  affect  the  result  would  be 
one  depending  on  the  cosine  of  the  declination,  or  of  the  form 

J/i'  —  k  cos  a. 

The  Pulkowa  catalogue  for  1865  is  the  best  modern  standard.  I  compared  it 
with  Auwebs's  declinations  for  1865,  on  which  his  proper  motions  rest,  with  a  view  of 
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detecting  any  error  in  the  latter  of  the  kind  in  question  The  discrepancy  between 
Auwers  and  Pulkowa  for  1865  for  each  octant  of  right  ascension  is  as  follows: 

Mean  differences  between  the  Pulkowa  declinations  for  1865  and  those  used  by  Auwers  for 

the  determination  of  proper  motion. 


Octant. 

P.-A. 
it 

Octant. 

P.-A. 
It 

I 

—  O.26 

V 

—  0.20 

II 

-O.28 

VI 

—  0.20 

III 

—  O.26 

VII 

-0.31 

IV 

-0-34 

VIII 

—  0.28 

It  is  clear  that  we  have  here  no  error  of  the  kind  sought  for.  If,  therefore,  there 
is  any  error  of  the  kind  sought  for  in  Auwers's  proper  motions,  it  must  arise  from 
some  other  cause  than  this. 

In  order  to  clear  up  the  matter  as  thoroughly  as  possible,  I  made  a  comparison  of 
Auwers's  proper  motions  with  my  own  new  standard  still  unpublished,  which  would  be 
free  from  any  hypothesis  whatever  as  to  the  way  in  which  Auwers  obtained  his  data. 
In  this  new  standard  catalogue  there  are  about  800  stars  whose  proper  motions  are 
found  in  Auwers'-Bkadley,  and  were  used  in  the  preceding  work.  I  made  a  compari- 
son of  these  proper  motions  with  those  of  Auwers.  The  difference  by  octants  is  as 
follows: 

Mean  corrections  of  Auwers's  proper  motions  as  given  by  the  standard  of  i8g6. 


Octant. 

N.-A. 

No.  of 

Octant. 

N.-A. 

No.  of 

11 

stars. 

II 

stars. 

I 

—  O.65 

114 

V 

—  0.20 

96 

II 

—  O.64 

I  12 

VI 

—  0.46 

IO9 

III 

—  O.7O 

98 

VII 

-0.53 

125 

IV 

—  O.78 

86 

VIII 

-0.47 

Il8 

It  is  evident  that  we  have  no  discordance  depending  on  the  right  ascension  which 
can  explain  the  discrepancy.  I  am  therefore  obliged  to  regard  the  discrepancy  in 
question  as  inexplicable,  and  to  reach  the  conclusion  that  no  other  result  than  that 
reached  in  the  present  paper  can  be  derived  from  the  Bradley  stars. 

Section  XVII. 

Precessionai  motions  derived  from  the  M.  A.'s  of  Lalande's  zones. 

The  recent  zone  catalogues  of  the  Astronomische  Gesellschaft  might  throw  much 
light  on  the  question  and  guide  us  in  adopting  a  mean  value  were  it  certain  that  the 
observations  on  which  they  are  based  were  comparable  with  the  ancient  ones.  In 
the  case  of  the  right  ascensions  there  is  little  difficulty  in  securing  this  comparability 
in  every  point  except  one,  the  personal  equation  dependent  on  magnitude.  The  prob- 
able effect  of  this  source  of  error  can  be  best  discussed  after  we  have  made  some 
comparisons. 

Of  the  older  material  the  portion  which  seems  best  worth  using  is  comprised  in 
the  zones  of  Lalande.    The  later  zones  of  Bessel  lose  value  not  only  on  account  of 
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their  more  recent  epoch,  but  from  the  fact  that  they  have  not  been  completely  com- 
pared with  a  modern  system.  In  the  case  of  two  zones  of  the  Gesellschaft  the  dis- 
cussion of  the  comparisons  with  Lalande  is  so  complete  that  no  difficulty  will  be 
found  in  studying  the  subject.  We  first  remark  that  the  Lalande  positions  as  now 
used  rest  on  Von  Astkn's  tobies,  which  again  depend  on  Piazzi's  positions  for  1800. 
It  is  therefore  to  be  presumed  that  the  systematic  corrections  to  reduce  Lalande's 
positions  to  a  modern  system  are  the  same  as  in  the  case  of  the  Piazzi  catalogue. 

Comparisons  of  the  older  catalogues  have  been  made  in  the  first  place  with  the 
A.  G  system.  There  is,  however,  a  certain  amount  of  doubt  as  to  definition  of  this 
system  as  adopted  by  different  investigators.  Since  the  system  in  question  is  com- 
monly supposed  to  be  based  on  a  comparison  of  the  Bradley  positions  as  determined 
by  Auwers  for  the  epoch  1755,  with  the  Pulkowa  observations  for  1865,  it  would 
naturally  be  supposed,  and  seems  to  have  been  very  generally  taken  for  granted,  that 
the  system  in  question  is  defined  by  this  comparison.  Such  is  not  the  case.  By 
carrying  back  the  A.  6.  positions  to  1755  we  do  not  reproduce  the  Bradley  places, 
but  a  system  of  places  in  the  general  average  about  o*.035  greater  in  R.  A.  This  arises 
from  the  fact  that  Auwers  derived  the  proper  motions  by  a  comparison  of  Bradley 
with  a  system  of  modem  positions  for  1865  derived  before  the  definitive  positions 
were  decided  upon.  When  a  decision  was  reached,  it  was  found  that  the  new  system 
gave  right  ascensions  about  o".035  greater  than  those  used  in  deriving  the  proper 
motions.  But  the  latter  were  used  without  change.  The  result  is  the  difference  in 
question.  Auwers  distinctly  states  this  fact  in  the  introduction  to  the  fundamental 
catalogue  in  Publication  XIV  of  the  A.stronomische  Gesellschaft.  But,  nevertheless, 
the  fact  seems  to  have  been  very  generally  overlooked.  The  result  of  this  state  of 
things  is  that,  in  the  absence  of  specific  statement,  we  must  be  in  doubt  as  to  which 
interpretation  has  been  employed.  I  shall,  however,  assume  that  Auwers  used  the 
A.  G.  system  and  not  the  Bradley  places. 

Beginning  now  with  Auwers's  Berlin  zone,  150  to  200,  we  find  on  pages  62-65 
a  very  exhaustive  comparison  with  Piazzi.  At  the  bottom  of  page  63  certain  correc- 
tions to  Piazzi's  right  ascensions  are  found,  which  make  them  comparable  with  the 
A.  G.  system.  The  values  vary  but  little  with  the  right  ascension,  and  the  mean 
throughout  the  circle  is 

^/  «  =  +  0*.202. 

This  we  take  as  the  correction  to  reduce  Lalande's  right  ascensions  to  the  A.  G.  sys- 
tem within  the  limits  of  the  zone  in  question 

Next,  on  pp.  230-246  of  his  catalogue,  Auwers  gives  a  comparison  of  Lalande's 
positions,  brought  forward  to  1875  by  precession  alone,  with  the  positions  observed 
by  him  in  his  zone.  I  have  taken  the  mean  of  the  excess  Beulin-Lalande  for  every 
hour  of  right  ascension,  omitting  stare  of  which  the  deviation  exceeded  one  second, 
and  also  those  marked  by  Auwers  with  an  asterisk  to  indicate  proper  motion.  The 
mean  value  of  the  excess  is 

J  a=r  +  o8.i75. 

This  agrees  well  with  the  mean  of  the  comparisons  given  by  Auwers  on  page  67, 
where,  however,  the  differences  are  given  by  zones  and  not  by  hours  of  right  ascension. 


RESULTS  OF  LALANDE'S  ZONES. 


59 


The  difference  of  the  two  corrections 


J  a  1&75  —  J  a  1800  =  — o".027 

is  to  be  taken  as  a  correction  to  the  precession  during  the  interval  between  the  obser- 
vations. This  differs  slightly  from  the  mean  of  the  sraoothed-off  results  given  by 
Auwebs  at  the  bottom  of  page  67,  which  would  be  alwut  —  o*.o6.  But  I  shall  use  it 
for  the  correction  to  Struve's  precession  in  74  years.  This  gives  a  correction  to 
Stbuve'8  ioow  of  —  o".$$.  But  the  centennial  motion  of  the  system  Nx  is  greater 
than  that  of  the  A.  G.  in  R.  A.  by  i".50.  Hence,  for  the  system  N,  we  should  1  ave 
as  the  correction  to  Struve's  m, 

100  J  m  =  +  °"-97- 

For  the  correction  to  the  provisional  values  forming  the  basis  of  the  present  inves- 
tigation, the  result  would  be, 

100  J  m  =  -\-  2". 27. 

This  is  an  entirely  inadmissible  result,  which  must  be  attributed  mainly  to  the  error 
depending  on  personal  equation  in  magnitude. 


Section  XVIII. 
Comparison  of  Boss's  Albany  zone  with  Lalande. 

In  the  Astronomical  Journal,  Volume  IX,  Boss  discusses  the  constants  of  solar 
motion  as  derived  from  the  comparison  of  his  zone  observations  (+ 1°  to  +5°)  with 
the  zones  of  Lalande,  Bessel,  etc.  Incidentally  he  refers  to  the  correction  of  the 
precessional  constant  which  may  thus  be  derived,  reaching  the  conclusion  that  correc- 
tions of  Struve's  m  and  »,  amounting  to  about  — o"-5  are  feebly  indicated. 

The  comparisons  with  Lalande  which  he  publishes  in  the  introduction  to  his 
zone,  pages  22-23,  lead  to  results  admitting  of  more  definite  statement.  He  gives  the 
following  summary  tobies  of  comparison  between  the  right  ascensions  of  Lalande, 
brought  forward  with  precession  alone  to  1875,  an(i  those  of  his  zone.  The  stars  are 
divided  into  two  classes,  bright  and  faint,  the  first  comprising  the  stars  up  to  the  sev- 
enth magnitude,  the  second  those  of  higher  magnitudes.  The  following  are  the 
results  in  each  sextant  of  the  circle  of  right  ascension: 


R.  A. 

o  h.  to  4  h. 
to  8 
to  12 
to  16 
to  20 
to  24 

Means 


4 
8 
1 2 
16 
20 


Mag 

Aa 
s 

.  1  to  7 
AS 

ff 

Mag. 
Aa 
s 

7.8  to  10. 
AS 

// 

Bright 
Aa 
s 

•  Faint. 
AS 

n 

+  O.I  I 

-3-6 

+  O.I4 

-3-8 

—  O.O4 

4-0.0 

+  O.I4 

-37 

4-O.I6 

-3-3 

—  O.O5 

—  0.2 

-O.08 

-2.7 

4-0.03 

—  2.2 

—  O.IO 

-O.4 

—  O.O3 

-37 

+  O.06 

-2.8 

—  O.O9 

—  O.I 

+  O.IO 

-2.7 

+  0.13 

-3  1 

—  O.O3 

-0-3 

+  O.I7 

-2-5 

4-0.20 

-4.0 

—  0.02 

+  O.9 

+  0.07 

-3-2 

4-O.I2 

-3-2 

—  O.05 

O.O 

6o 
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Iii  this  table  the  Albany  right  ascensions  are  not  corrected  for  personal  equation 
depending  on  magnitude.  When  thus  corrected  the  mean  Ja  for  the  group  i  to  7 
would  have  been 

•+•  o*.07  —  o".03  =  +  o'.04 ; 
and  the  mean  da  for  the  group  7.8  to  10  would  have  been 

-f  o".  1 2  —  o'.o7  =  +  o\<>5. 
The  mean  of  these  two  is  +  o".045. 

Next,  for  the  reduction  of  the  Lalande  positions  to  the  modern  system,  Boss 
formed  a  "special  catalogue,"  comprising  the  few  fundamental  stars  of  the  A.  G.  system 
which  lay  within  his  zone,  and  a  large  number  based  on  the  combination  of  Bradley's 
positions  with  the  modern  A.  G.  positions,  mistakenly  supposing  that  he  thus  has  a 
system  comparable  throughout  with  the  A.  G.  The  comparison  of  the  Lalande  cata- 
logue with  this  special  catalogue  is  as  follows  : 

R.  A.  Decl.  R.  A.  Decl. 


a 

Obs. 

s 

Obs. 

J8 

// 

a 

Obs. 

s 

Obs. 

AS 
tt 

oh. 

27 

+  O.I5 

22 

-3-2 

12k 

26 

+  O.I  I 

27 

-2.7 

2 

21 

+  O.IO 

18 

—  2.0 

14 

I  I 

+  O.I7 

I  I 

-2-4 

4 

H 

+  O.I6 

14 

-2.7 

16 

23 

+  0.20 

20 

-  1-5 

6 

22 

-j-O.21 

18 

—  2.  [ 

18 

13 

+  O.I6 

M 

-2.6 

8 

8 

+  O.OI 

8 

—  2.4 

20 

IO 

+  O.I9 

IO 

-3-9 

10 

16 

4-O.I9 

16 

-2-3 

22 

21 

+  O.II 

21 

-2.3 

The  general  mean  value  of  Ja  is  4-o"-I47-  Subtracting  this  from  o".045,  the 
modern  comparison  when  the  Albany  positions  are  corrected  for  personal  equation,  we 
have  a  difference  — o".io2  as  a  correction  to  the  general  precession  in  right  ascension 
during  a  period  of  eighty-four  years.  This  gives  as  the  correction  to  the  centennial 
precession  in  right  ascension  of  Struve  and  N0,  respectively,  the  amounts, 

s  " 

d  loom-  —  0.122  r: —  1.82  (corr.  to  Stbuve) 
4  roo  m  rr  —  0.034  —  —  0.52  (corr.  to  N„). 

But  the  general  mean  motion  in  right  ascension  of  the  system  Nx  is  greater  than 
that  of  the  Bradley-A.  G.  system  by  o".070.  We  have,  therefore,  relatively  to  the 
system  Nj  the  correction  of  the  general  precession  N0  in  right  ascension 

4  100  m—  1  ".05  —  o".52  =  +  °"-53- 

We  may  accept  this  as  the  definitive  correction  given  by  the  comparison  of  Boss 
and  Lalande. 

I  have  made  no  attempt  at  an  independent  investigation  of  the  comparisons  in 
declination,  because  of  the  unknown  possible  systematic  error  of  the  Lalande-Piazzi 
declinations.  It  may  be  remarked,  however,  that  the  correction  of  -o"-5  given  by  Boss 
to  Struve's  precession  in  declination  is  eauivalent  to  a  correction  of  the  provisonal 
value  of  ioow,  amounting  to 

// 

J  100  »  =  -1-0.34. 
J  100  w  =  + 0.80. 
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The  interesting  feature  of  these  results  is  that  they  bring  the  correction  of  m  more 
nearly  into  accordance  with  n.  But  in  accepting  the  results  from  the  right  ascension 
as  real,  we  meet  the  difficult  question  of  the  influence  of  magnitude  on  personal 
equation.  My  investigations  on  this  subject,  which  will  appear  with  some  fullness 
in  the  next  part  of  the  Astronomical  Papers,  and  of  which  a  summary  is  given  in  the 
Astronomical  Journal,  Volume  XVI,  p.  65,  lead  to  the  conclusion  that  in  the  case  of 
the  leading  principal  observatories  this  personal  equation  has  nearly  the  same  value 
for  eye  and  ear  observations  that  it  has  for  those  registered  011  the  chronograph. 
Although  this  result  was  unexpected,  and  is  probably  contrary  to  the  general  impres- 
sion, we  have  a  very  simple  explanation,  nearly  identical  with  that  now  accepted  as 
the  cause  of  the  corresponding  error  in  the  case  of  the  chronographic  observations. 
In  both  cases  the  observer  is  biased  in  his.  estimate  of  the  distance  of  the  star  from 
the  transit  thread  by  the  apparent  magnitude  of  the  star.  When  the  latter  is  apparently 
smaller,  his  unit  of  measurement  is  unconsciously  taken  less  in  proportion,  and  the 
distance  from  the  wire  is  therefore  assumed  to  be  greater.  In  the  case  of  eye  and  ear 
observations  this  erroneous  estimate  before  the  star  crosses  the  thread,  should,  in  the 
general  mean,  be  neutralized  by  the  corresponding  error  on  the  opposite  side  after 
the  star  has  crossed  the  thread.  Were  the  position  of  the  star  at  the  first  clock-beat 
following  transit  estimated  by  the  observer  with  the  same  leisurely  precision  as  at  the 
moment  of  the  clock-beat  next  preceding,  it  is  to  be  supposed  that  the  error  would  not 
exist.  But,  in  the  case  of  modern  observations,  when  the  observer  only  has  fifteen  or 
twenty  seconds  between  the  several  threads  in  which  to  record  his  observation  and  pre- 
pare for  the  transit  over  the  thread  next  following,  he  is  naturally  hurried  and  fails  to 
make  a  sufficiently  careful  estimate  after  the  star  has  passed  the  thread.  His  result  is 
therefore  affected  by  the  bias  of  the  approaching  star. 

In  the  case  of  the  older  transits,  the  observations  were  either  made  over  a  single 
thread,  or  the  threads  were  so  wide  apart  that  the  observer  had  plenty  of  time  for  his 
estimate.  It  is  therefore  quite  likely  that  in  these  cases  the  error  in  question  does  not 
exist  A  comparison  of  Boss-Lalande  seems  to  show  that  it  did  not  exist  in 
Lalande's  observations.  We  must,  not,  however,  too  hastily  conclude  that  by  the 
simple  device  of  observing  a  star  over  a  single  thread,  or  by  taking  plenty  of  time  for 
the  estimate,  the  error  would  necessarily  be  avoided.  There  might  be  fear  that,  in 
any  case,  the  estimate  made  after  the  star  had  crossed  the  thread  -would  be  biased  by 
the  preceding  expectation  of  the  observer.  This  bias  ought,  however,  to  be  done 
away  with  through  proper  discipline  on  the  observer's  part. 

Returning  to  the  two  preceding  values,  since  the  error  in  question  was  absolutely 
determined  by  Boss  in  the  case  of  his  zone  and  is  therefore  eliminated,  we  may  regard 
the  result  of  his  comparison  with  Piazzi  .as  one  not  affected  by  this  error.  In  the  case 
of  Schjellebup  the  presumption  would  be  in  favor  of  the  existence  of  the  error. 
The  result  seems  to  show,  however,  that  it  was  very  small.  But  we  can  not,  on  this 
ground  alone,  assume  that  it  vanished,  and  therefore  assign  a  large  weight  to  his 
result. 
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Section  XX. 

Elimination  of  the  parallactic  motion  from  the  precession  of  each  individual  star. 

There  is  yet  a  fourth  mode  of  treating  the  subject,  which  is  eminently  worthy  of 
consideration,  although  I  am  unable  to  apply  it  exhaustively.  If  we  regard  the  posi- 
tion of  the  solar  apex  as  known,  then  we  may  eliminate  the  parallactic  motion  from  the 
two  equations  given  by  each  particular  star,  and  thus  obtain  an  equation  in  which  the 
precession  appears  as  the  only  unknown  quantity.  This  amounts  to  the  same  thing 
as  determining  the  magnitude  of  the  precession  by  the  condition  that  the  apparent 
rotation  of  the  whole  system  of  stars  around  the  direction  of  the  solar  apex  as  an  axis 
shall  vanish.  The  apparent  motion  of  precession  may  be  considered  as  one  taking 
place  around  the  pole  of  the  ecliptic  as  an  axis.  It  is  therefore  a  fortunate  circum- 
stance in  the  application  of  this  method  that  the  distance  between  the  solar  apex  and 
the  pole  of  the  ecliptic  is  only  about  300.  The  result  is  that  our  equations  do  not 
lose  much  in  weight  by  tho  elimination  in  question.  The  latter  is  conducted  in  the 
following  way.    Let  us  put 

A,  D ;  the  R  A.  and  Dec.  of  the  solar  apex, 
q ;  the  parallactic  motion  of  a  star. 

We  then  have  for  the  components  of  the  solar  motion, 
X  —  pq  cos  D  cos  A 

Y  —  pq  cos  D  sin  A 
Z  =  pq  sin  D 

We  also  put 

a'  =  cos  D  sin  (a  —  A) 

V  =  cos  D  sin  6  cos  (a  —  A)  —  sin  D  cos  S 

By  substitution  of  these  values  of  X,  Y,  and  Z,  in  the  equations  of  p.  33  they  become 

aSp  -\-dq  —  n  cos  8 
Up  +  b'q  =  m' 

Eliminating  q  from  these  equations  and  putting  for  brevity 

n-a'b  —  ab' 

we  have 

IISp  —  a'  n' —  b'/i  cos  8 

By  substituting  in  IT  the  values  of  a,  a',  b,  and  V,  we  see  that  it  may  be  thrown  into 
the  form 

A  sin2  6  +  B  sin  S  cos  6  +  C  cos2  S  +  E 
and  this  again  into  the  form 

L  +  M  sin  2  S  -\-  N  cos  2  6 

the  coefficients  L,  M,  and  N  being  functions  of  a.  Thus,  a',  V ,  and  IT  may  be  tabu- 
lated as  functions  of  a  and  6* 

•These  formula;  are  not  the  best  for  the  purpose,  since  they  ignore  the  fact  that  near  the  solar  apex  the  paral- 
lactic motion  vanishes.  As  the  result  derived  in  this  section  does  not  enter  directly  into  the  conclusion,  I  deem  it 
unnecessary  to  study  the  best  formulae  at  present 


64 


THE  PRECESSIONAL  CONSTANT. 


Our  problem  requires  the  determination  of  the  values  of  A  and  D,  which  dete- 
mine  the  position  of  the  solar  apex.  The  determination  of  this  quantity  which  we 
have  already  made  should  not  be  regarded  as  the  definitive  one,  because  it  depends 
solely  on  stars  with  small  proper  motions.  If  the  stars  had  no  absolute  proper  motion, 
it  is  clear  enough  that  the  solar  motion  could  best  be  determined  from  the  stars  having 
the  largest  parallactic  motion,  because  they  would  be  the  nearer.  In  any  case,  in  the 
general  mean,  the  presumption  is  that  the  stars  with  a  larger  proper  motion  are  the 
nearer  to  us  and  therefore  best  adapted  to  the  determination  in  question.  At  the  same 
time  the  method  which  has  always  been  pursued  of  selecting  stars  having  large  proper 
motions  and  determining  the  parallactic  motion  from  them  alone  is  not  free  from  objec- 
tion. The  determination  of  the  proper  motion  of  a  star  presupposes  a  knowledge 
of  the  magnitude  of  the  precession;  and  a  star  may  or  may  not  be  included  in  the 
list  of  those  having  a  sufficient  proper  motion,  according  to  whether  one  value  or 
another  of  the  precession  is  used. 

It  seems  to  me  that  the  numerous  results  obtained  by  past  investigators  should 
now  give  way  to  the  three  latest  determinations  by  L.  Struve,  Bosh,  and  Stumpe. 
From  the  Bradley  stars  L.  Struve  found: 

A  =  2730.4;  D  =  270.3. 

But,  owing  to  the  systematic  correction  to  Bradlev's  declinations  which  I  have  found, 
I  estimate  that  his  parallactic  motion  Z  requires  a  correction  of  -\-o".yo  to  reduce  it  to 
the  standard  which  I  have  called  C^.  With  this  correction  we  may  put  the  results 
into  the  form 

X  =  +0.23 

Y  =  -3-87 
Z  =  +  2.70 

This  gives,  as  L.  Struve's  corrected  result, 

o 

A  =  2734 
D=  34-9 

Boss's  general  result  (A.  J.,  Volume  IX,  page  1 65)  is 

o 

A  —  280 
D  -  40 

Stumpe,  in  his  very  exhaustive  discussion  (Astronomische  Nachrichten,  volume 
140,  pages  177-190),  gives  a  number  of  separate  results  which  have  the  great  advan- 
tage of  being  derived  from  southern  stars  as  well  as  northern  stars.  I  have  made  a 
general  combination  of  them  to  derive  what  seems  to  me  the  most  probable  final  result 
from  his  work,  and  find  it  to  be, 

o 

A  =  279.5 
D=  38.7 
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Altogether  it  seems  to  me  that  we  may  take  the  most  likely  position  of  the  solar 
apex  as, 

o 

A  =   277.5  =  i8k  30m 
D  =  +38.0 

This  conclusion  places  the  solar  apex  in  the  middle  of  the  eighteenth  hour  of 
right  ascension,  in  a  direction  different  from  that  of  the  star  a  Lyras  by  a  quantity 
not  exceeding  its  probable  error. 

After  preparing  tables  based  on  this  adopted  position  of  the  Rolar  apex  and  on 
the  preceding  formula?,  I  find  that  these  formula;  are  not  those  best  adapted  to  the 
work,  and  that  they  coincide  with  the  best  formula?  only  near  what  we  may  call  the 
sun-way  equator,  or  the  belt  900  distant  from  the  solar  apex.  In  all  other  positions, 
the  measures  being  considered  as  those  of  rotation  around  the  solar  apex,  it  will  follow 
that  each  star  enters  into  the  result  with  a  weight  proportional  to  the  cosine  of  the 
sun-way  latitude,  or  to  the  sine  of  its  distance  from  the  apex.  At  the  apex  itself  the 
rotation  indicated  by  the  proper  motion  necessarily  enters  the  equation  with  a  van- 
ishing weight,  and  therefore  gives  no  result  whatever.  But  it  is  evident  that,  since 
the  parallactic  motion  vanishes  at  the  solar  apex,  all  stars  near  this  apex  may  be 
included  in  the  result  without  the  necessity  for  eliminating  the  parallactic  motion.  It 
is  true  that  since  such  stars  are  only  about  300  distant  from  the  pole  of  the  ecliptic 
they  will,  in  any  case,  only  enter  with  a  small  weight,  yet  they  should,  in  theory,  be 
included.  The  method  in  question  is  substantially  equivalent  to  the  determination  of 
the  precession  from  the  proper  motion  of  the  star  in  sun-way  longitude  only.  I  find 
that  this  method  of  determination  was  proposed  several  years  ago  by  Kapteyn,  but  I 
became  aware  of  it  too  late  to  utilize  the  idea  in  the  form  in  which  he  presents  it. 
It  differs  from  my  form  only  in  the  practical  introduction  of  a  different  system  of 
weights  for  stars  distant  from  the  sun-way  equator. 

My  method  of  proceeding  was  determined  by  the  following  considerations.  The 
parallactic  motion  is  not  only  eliminated  from  the  precession  in  the  immediate  neigh- 
borhood of  the  apex,  where  it  vanishes,  but  it  is  eliminated  from  the  proper  motion  in 
right,  ascension  on  both  the  meridian  and  the  antimeridian  of  the  apex.  That  is  to 
say,  from  proper  motions  in  right  ascension  of  stars  near  these  meridians  we  have 
a  determination  of  the  precession  from  which  the  solar  motion  is  eliminated. 

I  consider  that  all  stars  within  one  hour  of  the  apical  meridian  and  antimeridian 
satisfy  this  condition,  for,  although  at  a  distance  of  150  from  that  meridian,  the  effect 
of  the  parallactic  motion  will  be  quite  appreciable ;  it  would  be  in  opposite  directions 
on  the  opposite  sides  of  the  meridian,  and  therefore  ought  to  be  eliminated  from  the 
mean  result.  Taking  out  the  two  lunes  thus  defined  from  the  sphere,  we  have  left 
two  regions,  the  one  extending  in  R.  A.  from  I9b.5  through  the  vernal  equinox  to  5h.5; 
the  other  from  7b.5  through  the  autumnal  equinox  to  17*. 5.  I  shall  call  these  two 
regions  the  vernal  and  autumnal  regions,  respectively.  Within  them  I  took  all  the 
stars  which  were  within  300  of  the  sun-way  equator  or,  more  exactly,  all  those  for 
which  the  value  of  the  coefficient  IT  exceeded  0.80,  and  from  the  proper  motions  in 
each  region  a  value  of  the  precession  was  derived  by  the  formula;  which  should  be 
independent  of  the  parallactic  motion. 
A  p,  VOL  VIII,  pt  1  5 
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In  the  case  of  the  apical  lunes  I  simply  took  the  normal  equations  in  right  ascen- 
sion as  already  formed  for  all  the  individual  stars,  except  those  of  the  first  and  second 
magnitudes.  A  combined  equation  for  each  lune  was  then  formed  from  the  separate 
ones  by  the  factors  for  magnitudes  already  given  with  the  following  results: 

Apical  lune;        1558  Sp  —  162 iX  —  258Y  =  +  1826" 
Antiapical  lune ;    2940  Sp  +  2674X  +  2 1 6Y  =  4-  2367 

From  each  of  these  combined  equations  we  may  obtain  a  result  in  two  ways.  In 
one  way  we  assume  that  the  parallactic  motion  is  eliminated  from  the  mean,  and 
simply  ignore  X  and  Y.  In  the  other  way  we  substitute  the  values  of  X  and  Y  as 
already  found  for  the  standard  of  the  fifth  magnitude  stars.  Beginning  with  the  apical 
lune,  when  we  omit  X  and  Y,  we  have 

^=1826:  1558  =  +  1".  1 7. 

Retaining  X  and  Y,  the  absolute  term  becomes 

1826  +  535  —  650=  171 1, 

whence 

<Jp  =  +i".io. 

Passing  to  the  antiapical  meridian  we  have,  when  we  omit  X  and  Y, 

<Jp  =  4-o".8i. 

Retaining  them,  the  second  member  of  the  equation  becomes 

2367  —  543  +  544  =  2368, 

whence  we  have  the  same  value,  +o".8i. 

The  discordance  of  the  results  on  the  apical  and  the  antiapical  lunes  would  be 
removed  by  taking  a  smaller  right  ascension  of  the  solar  apex,  say,  the  value  already 
found  in  this  investigation.  But  I  do  not  consider  that  this  affords  any  sound  evidence 
for  doubting  the  superiority  of  the  adopted  value.  Table  XI  shows  the  computation 
of  the  motions  in  the  vernal  and  autumnal  regions,  bounded  by  the  apical  lunes. 


Table  XI. — Computation  of  the  precaution  from  902  start  near  the  ecliptic  and  the  gun-way  equator 
by  eliminating  the  parallactic  motion  from  the  equation*  given  by  each  ttar. 


Vernal  region;  I^.s  to  5h.5  R.  A. 

Autumnal  region ;  7h 

•5  to  I7h.s  R.  A. 

Mag. 

Stars. 

Sr  sin  J 

Prod,  by  KJ. 

Sp. 

. 

Mag.  Stars.  2rs\nd 

1  ! 

211 

Prod,  by  F*. 

«P 

// 

n 

// 

// 

// 

3 
4 
5 
6 

7 

18 
64 
149 
231 
56 

+  15- 4 
46.0 

56.5 
189.3 
56.2 

15-  5 
55-4 
128.3 

195-7 
48.6 

5-5  5-5 
29-4  35-5 
56.5  128.3 
272. 6  281.  8 
1 10.  2   95. 3 

+°-99  : 

+0.83 . 

+0.44 
40.97 
+1. 16  1 

3 
4 

5 
5 
7 

19 
44 
120 

214 
47 

-'3-7 
+  a  1 
+24.0 
+  3-9 
-  3-4 

16.  1 

37-7 
100. 1 
182.  2 

4°- 5 

—  4-9  5-8 
+  0.  1    24.  r 
+24. 0  100.  1 
+  5-6  58.3 

-  6-7  79-4 

—0.85 
0. 00 
+0.24 
+0. 02 
—0.08 

5i8 

474-  2  546.  4 

+0. 87 . 

444 

18. 3  267.  7 

-fo.07 
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It  will  be  seen  that  the  stars  have  been  first  classified  according  to  the  magnitudes, 
as  in  the  former  portions  of  the  work.  In  the  third  column  the  symbol  r  may  be  con- 
sidered as  representing  the  proper  motion  of  the  star  on  a  great  circle  resolved  in  the 
direction  of  sun-way  longitude,  d  represents  the  distance  of  the  star  from  the  solar 
apex.  The  combination  t  sin  d  is  then  the  product  of  the  proper  motion  by  the 
appropriate  factor  when  a  solution  is  made  by  least  squares  in  order  to  represent  r  sin  -J 
by  the  processional  motion. 

The  next  column,  2TI,  contains  the  sum  of  the  corresponding  coefficients  of  the 
lunisolar  precession  when  normal  equations  for  this  precession  are  formed. 

In  order  to  form  the  final  normal  equation  these  partial  normals  are  multiplied  by 
the  squares  of  the  factor  already  found,  which  are  taken  to  represent  the  relative  mean 
distances  of  the  stars  of  the  several  magnitudes.    The  factors  then  are,  respectively, 

0.6,  0.8,  1.0,  1.2,  and  1.4. 

There  is  an  evident  systematic  difference  in  the  results  from  the  two  regions,  the 
largest  result  from  the  autumnal  region  being  smaller  than  the  smallest  result  from 
the  vernal  region.  This  difference  extends  through  all  the  magnitudes,  and  can  there- 
fore not  be  considered  as  accidental. 

For  the  sake  of  study  we  place  the  four  results  just  derived  from  stars  in  corre- 
sponding but  not  identical  regions,  in  juxtaposition  with  those  already  derived  from 
the  right  ascensions  alone.  In  the  case  of  the  latter  I  have  taken  the  results  already 
given  for  the  octants  including  the  apex  and  the  antiapex,  and  the  weighted  means  of 
the  intermediate  octants: 

From  motion  in 
From  R.  A.  sun-way  longitude. 

//  // 

Apical  region  .  .  .  dp  =  +1.04  -f-1.14;  wt  =  i 
Vernal  region  .  .  .  +0.30  -(-0.87;  2 
Antiapical  region  .  .  +  0.40  4-°«8i;  1 
Autumnal  region  .    .         —0.21       +0.07;  2 

Giving  double  weight  to  the  results  from  the  vernal  and  autumnal  regions,  the 
number  of  stars  being  nearly  double,  the  result  of  eliminating  the  parallactic  motion  is 

«Jp=+o".64. 

As  to  the  R.  A's  I  have  already  indicated  that  the  discordances,  or  the  small 
values  given  by  the  intermediate  octants,  may  be  accounted  for  by  supposing  a  system- 
atic or  accidental  difference  of  absolute  magnitudes  in  opposite  quarters  of  the  sphere. 
Were  such  the  case,  any  such  difference  should  be  eliminated  with  the  solar  motion,  and 
therefore  not  appear  in  the  results  derived  after  that  motion  is  eliminated.  It  will  be 
seen  that  the  discordance  is  removed  in  the  case  of  the  vernal  regions,  but  not  in  that 
of  the  autumnal  regions.  It  is  also  remarked  that,  supposing  this  difference  of  magni- 
tude to  be  real,  it  affords  an  explanation  of  the  larger  value  found  for  the  declinations. 
A  very  little  consideration  will  show  that  if  the  stars  of  a  given  apparent  magnitude 
are  farther  away  within  the  vernal  region  than  within  the  autumnal  region,  then  the 
smaller  parallactic  motions  in  the  former  region  will  tend  to  diminish  the  precession 
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found  from  the  right  ascensions,  and  increase  that  found  from  the  declinations.  The 
same  effect  will  be  produced  by  the  relatively  greater  negative  parallactic  motions  of 
the  nearer  stars  in  the  autumal  region. 

But  in  view  of  the  almost  uniform  extension  of  the  discordance  through  all  the 
apparent  magnitudes,  as  shown  by  a  comparison  of  the  colums  Sp  in  Table  XI,  I 
still  think  the  main  cause  of  discordance  to  be  an  error  of  double  period  in  the 
Bradley  R.  A.'s. 

Section  XXI. 

Summary  and  discussion  of  results. 

The  various  results  for  the  precessional  motion  derived  and  discussed  in  the  pre- 
ceding pages  rest  on  two  independent  bases,  the  observed  motions  in  right  ascension, 
and  those  in  declination.  The  motion  in  sun-way  longitude  does  not  rest  on  a  third 
independent  basis,  but  should  rather  be  considered  as  a  particular  combination  of  the 
other  two  motions.  What  we  have  finally  to  do  is  to  decide  independently  upon  the 
most  probable  result  of  the  motions  in  right  ascension  alone,  then  on  the  most  probable 
result  from  the  declinations  alone,  and  then  upon  the  best  combination  of  the  two. 

Taking  first  the  right  ascensions,  we  have  found  from  the  Bradley  stars,  by 
different  modes  of  treatment,  the  following  values  for  correction  to  the  value  N„  of  the 
lunisolar  precession : 

// 

I.  From  individual  stars,  each  with  equal  weight  «$p  =  -J-o.24 

II.  From  groups  of  stars,  each  having  an  adjusted  weight;  +0.20 

III.  By  division  into  lunes,  each  an  octant  in  extent,  giving 

greater  weight  to  the  lunes  nearer  the  meridian  of  the 

solar  apex ;  +  o  30 

IV.  By  the  statistical  method;  +°-37 

Of  these  four  results  I  conceive  that  the  third  is  best  entitled  to  acceptance  as  the 
definitive  value  to  be  derived  from  the  motion  of  the  Bradley  stars  in  right  ascension. 
While  the  results  of  the  three  other  methods  are  not  to  be  left  wholly  out  of  consider- 
ation we  can  not  say  that  their  combined  result  would  be  decidedly  different  from  o".30, 


some  being  greater  and  some  less. 

From  the  Lalande  stars  we  have  found  : 

By  comparison  with  Boss's  Albany  zone ;   +0-57 

By  comparison  with  Schjellerup;   +067 

Mean,   +0.62 


Of  the  two  results  here  combined  I  conceive  that  there  can  be  no  doubt  of  the 
legitimacy  of  including  Boss's  result,  since  the  injurious  effect  of  variation  of  personal 
equation  with  magnitude  was  carefully  determined  and  eliminated  from  the  compari- 
son of  his  right  ascensions  with  those  of  Lalande.  The  one  weak  point  is  that  his 
zone  was  only  40  wide.  Against  this  is  to  be  placed  the  consideration  that  the  large 
number  of  stars  observed  by  Lalande  within  this  zone  were  most  carefully  reduced 
by  comparison  with  a  sufficient  number  of  standard  stars  in  the  zone  itself.  On  the 
other  hand  we  have  no  evidence  that  Sch.iellerup's  right  ascensions  are  free  from  this 
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Were  it  possible  to  assign  any  law  of  the  planetary  precession  for  the  past  and 
future  in  this  way,  this  proceeding  might  be  worthy  of  consideration.  But  in  view  of 
the  necessary  uncertainty  which  attaches  to  the  homogeneity  of  the  system  Nx  itself, 
and  of  the  impossibility  of  forming  any  hypothesis  of  the  past  and  future  value  of 
the  planetary  precession  otherwise  than  on  the  basis  of  theory;  considering  also  the 
undoubted  character  of  the  theory,  and  the  fact  that  no  admissible  changes  in  the 
masses  of  the  planets  can  change  the  amount  of  this  motion  by  more  than  one  or  two 
tenths  of  a  second  per  century,  I  conceive  that  this  mode  of  proceeding  is  inadmissible. 
We  should  rather  anticipate  that  the  astronomy  of  the  not  distant  future  will  reconcile 
the  two  motions  with  theory. 

The  second  method  is  to  consider  J~EJ  as  an  unknown  correction  of  which  the 
most  probable  value  is  zero,  but  which  has  a  probable  error,  to  be  determined  as  best 
we  can.  This  probable  error  should  be  combined  with  that  of  the  absolute  term 
o".36.  The  question  would  then  merge  itself  into  that  of  the  probable  error  to  be 
assigned  to  the  value  of  p  as  derived  from  the  right  ascensions  alone.  Independently 
of  ^/E'  we  have  found  that  the  values  of  p  derived  from  pairs  of  lunes  of  opposite 
right  ascensions,  from  the  mean  of  which  the  parallactic  motion  is  eliminated,  differ 
by  more  than  i".  This  would  seem  to  imply  a  probable  error  of  not  less  than  o".$o 
in  the  absolute  term  of  p.  But  if  we  regard  the  point  of  reference  from  which  p  is  to 
be  measured  as  the  actual  mean  of  all  the  stars  around  the  circle  of  right  ascension, 
this  discrepancy  becomes  simply  a  periodic  term  in  the  correction  of  the  proper 
motions,  and  no  longer  indicates  a  probable  error  of  the  absolute  term. 

Passing  now  to  ^E',  under  the  hypothesis  at  present  under  consideration  its 
probable  error  would  be  simply  its  probable  value,  which,  in  the  absence  of  any 
evidence  as  to  that  value,  we  might  estimate  at  ±  o".20.  On  the  whole  we  may  con- 
sider that  the  range  of  judgment  as  to  what  probable  error  should  be  assigned  to  the 
value  of  p  as  derived  from  the  right  ascensions  is  quite  large,  and  that  it  might  be 
placed  anywhere  between  ±o".20  and  ±o".40. 

In  the  case  of  the  value  derived  from  the  declinations  there  are  fewer  sources  of 
uncertainty.  It  may  be  placed  between  ±  o".$o  and  ±  o"-40.  Considering  the  com- 
bination, there  seems  to  be  a  possible  range  of  doubt  whether  wo  should  regard  the 
result  from  the  declinations  as  entitled  to  more  weight  than  the  right  ascensions,  or  as 
entitled  only  to.  one-third  the  weight.  Altogether  it  would  seem  that,  from  this  point 
of  view,  we  can  only  say  that,  according  to  different  judgments,  the  value  of  Sp  might 
range  between  limits  -t-o".50  and  +o".8o,  or 

+  o".5o<dp<  +  o".8o. 

The  third  method  of  procedure  is  to  anticipate  the  probable  correction  which 
must  hereafter  be  applied  to  the  motion  of  the  equinox,  or  the  probable  value  of  JYJ. 
Although  a  definitive  determination  of  this  correction  is  not  possible  at  the  present 
moment,  there  are  data  for  an  approximate  estimate.  The  Greenwich  observations  of 
the  Sun  during  the  sixty  years  since  1 834  are  of  great  value  for  this  purpose,  owing 
to  their  having  been  made  on  a  uniform  system  by  so  large  a  number  of  observers 
that  the  personal  equation  peculiar  to  the  limb  of  the  Sun  is  probably  nearly  elimi- 
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nated  from  the  final  result.  I  find  that  taking  the  Sun's  absolute  longitude  as  given  in 
my  new  theory,  a  rough  investigation  of  the  Greenwich  results  during  the  sixty  years, 
1835-1895,  leads  to  the  approximate  value, 

JE'  =  +  o".5. 

I  also  conceive  that  the  observed  right  ascensions  of  Mercury  are  worthy  of  con- 
sideration in  this  respect,  because  they  are  nearly  free  from  the  effect  of  the  personal 
error  affecting  the  limbs  of  the  Sun.  It  is  true  that  the  Conference  of  1 896  expressed 
a  preference  for  determining  the  equinox  from  the  Sun  alone,  owing  to  the  manner  in 
which  uncertain  elements  of  Mercury  would  enter  into  a  result  derived  from  observa- 
tions of  that  planet.  I  conceive  that  this  objection  is  of  force  only  as  regards  the 
determination  of  the  Sun's  absolute  longitude  from  the  declination  of  Mercury. 
When  observations  of  right  ascension  alone  are  considered,  the  reduction  of  an 
observed  right  ascension  of  Mercury  to  the  center  of  its  motion  is  one  all  the  errors  of 
which  will  be  nearly  eliminated  in  any  one  revolution  of  Mercury  in  its  relative  orbit. 
I  conceive  therefore  that  the  correction  JE'  =  i".o  given  by  observations  of 
Mercury  alone  is  one  legitimately  worthy  of  consideration. 

These  corrections  are  not  to  be  accepted  by  themselves  alone,  but  are  to  be 
considered  in  connection  with  the  great  mass  of  material  on  which  the  equinox  for  the 
Bystem  Nj  is  based.  This  material  of  course  gives  the  value  zero  for  the  correction 
in  question.  Altogether  I  conceive  that  the  existing  material  justifies  us  in  anticipat- 
ing for  this  correction  a  value 

JE'  =  +  o"  30. 

Accepting  this  result,  we  shall  have  from  the  right  ascensions  alone, 

Sp  —  o".$6  -f- 1 .09  JE'  =  +  o".6g. 

The  combination  with  the  declinations  will  then  lead  to  a  probable  result  somewhere 
between  the  limits 

o".8o<($p<o".90. 

We  have  thus  a  correction  +o".8o  which  appears  to  be  about  the  upper  limit  of 
probability  if  we  regard  the  system  N,  as  that  to  which  our  equinox  is  to  be  referred, 
and  the  lower  limit  in  the  event  of  the  correction  which  now  seems  most  probable 
having  to  be  applied  in  the  future. 

In  what  precedes  we  have  considered  the  question  of  the  relative  weights  to  be 
assigned  to  the  results  from  the  right  ascensions  and  declinations  as  the  result  of  a 
somewhat  uncertain  judgment.  There  is,  however,  one  principle  which  will  lead  us 
to  a  fairly  definite  combination,  although  I  have  not  investigated  it  rigorously.  I 
have  already  remarked  that  a  possible  cause  for  the  discrepancy  between  the  two 
results  is  to  be  found  in  an  inequality  of  the  average  magnitude  of  stars  in  the 
autumnal  and  vernal  regions  of  the  heavens,  and  a  consequent  inequality  in  the  paral- 
lactic motion  of  stars  of  equal  apparent  magnitude  in  the  two  directions.  This  cause 
will  produce  opposite  effects  on  the  precessional  motions  as  derived  from  right  ascen- 
sion and  from  declination.  The  amount  of  the  effect  does  not  admit  of  absolute  state- 
ment, because  it  would  depend  upon  the  distribution  of  the  inequality  in  different 
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parts  of  the  heavens.  Now,  the  elimination  of  the  parallactic  motion  from  the  equa- 
tions given  by  each  star  eliminates  this  error  completely,  and  at  the  same  time  affords  a 
criterian  for  the  relative  weights  of  the  R.  A.'s  and  the  Dec's.    We  have  found  in  this  way 

fy  =  +  o".64 

as  that  of  the  combination  in  question. 

A  correction  of  +  o".30  to  all  the  motions  in  right  ascension  would  increase  this 
result  by  a  quantity  which  I  have  not  accurately  determined,  but  which  would  not 
differ  much  from  o".20.  With  this  increment  we  should  therefore  have,  by  eliminating 
the  parallactic  motion — 

Sp  =  +  o".84. 

The  question  now  arises,  whether  it  is  advisable  to  anticipate  this  correction  to 
the  motion  of  the  equinox.  It  seems  to  me  that  it  is,  for  these  reasons:  It  is  to  be 
expected  that  the  system  of  right  ascensions  will  be  corrected  from  time  to  time  to 
bring  it  into  accordance  with  the  observed  equinox.  Such  corrections  will  be  product- 
ive of  little  trouble  or  confusion.  But  a  change  in  the  precessional  motion  will  be  pro- 
ductive of  confusion  to  astronomers,  while  a  small  error  in  its  amount  will  not  cause 
any  trouble.  It  may,  therefore,  be  anticipated  that  such  value  of  the  precessional 
motion  as  may  now  be  generally  adopted  will  be  continued  through  a  largo  part  of 
the  twentieth  century,  perhaps  the  whole  of  it,  should  no  large  correction  be  found 
necessary. 

From  all  these  considerations  I  have  considered  that  the  correction 

Sp  =  +  o".8o 

would  be  very  near  the  best  compromise  to  adopt.  This  is  very  nearly  o.oooi  6  of  the 
lunisolar  precession  itself.  I  have,  therefore,  increased  the  provisional  value  of  the 
processional  constant  by  this  fraction  of  its  whole  amount,  and  reconstructed  the  pre- 
cessional motions  accordingly.  I  find  that,  owing  to  the  defects  of  decimals  of  the 
computations  employed  in  "Elements  and  Constants,"  the  actual  correction  thus 
derived  to  the  lunisolar  precession  there  given  is 

Sp  =  +  o".82. 

I  therefore  propose  this  correction  as  the  value  which,  considering  all  the  data  for 
judging  the  future,  seems  most  likely  to  meet  the  requirements  of  the  case.  If  it  be 
too  small,  the  injurious  effect  on  n  will  apparently  only  be  about  o".  1 2  per  century, 
even  if  we  assign  to  the  latter  the  value  derived  from  the  decimations  alone.  If  it 
prove  to  be  too  large,  the  excess  will  probably  not  produce  any  confusion  in  astronomy, 
since  it  will,  in  any  case,  be  but  a  small  fraction  of  the  parallactic  motion.  We  therefore 
have  as  the  definitive  conclusion  of  the  preceding  investigation  the  following  value  of 

the  precessional  constant  as  I  have  defined  it  and  of  the  precessional  motions  for  1850: 

// 

.    Precessional  constant;  P— 5490.66  —  o".oo364T 
.    General  precession;  ^=5024.53 
Lunisolar  precession;   ^  =  5036.84 
100  in  —  4607. 1 1 
100  n—  2005.1 1 
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With  this  value  of  the  precessional  constant  I  have  gone  over  the  investigation 
of  the  motion  of  the  equator  found  on  pages  196-202  of  my  "Elements  and  Con- 
stants," with  some  variations  of  methods,  designed  to  make  the  whole  more  rigorous 
and  homogeneous.    The  results  are  given  in  the  next  section. 


Section  XXII. 


Numerical  valves  of  the  precessional  motions. 

The  following  tables  show  the  numbers  to  be  used  in  reductions  of  the  mean 
places  of  the  fixed  stars  from  one  epoch  to  another  when  we  use  the  preceding  value 
of  the  precessional  constant,  and  the  secular  motion  of  the  ecliptic  as  deduced  from  the 
masses  of  the  planets  employed  in  my  tables  of  the  Sun.  The  data  are  of  two  classes, 
those  which  express  the  annual  motions  of  precession,  and  those  serving  for  the  rigor- 
ous trigonometric  reduction  of  star  places  from  one  epoch,  to  another. 


Values  of  the  centennial  precessional  motions  from  1725  to  2000. 

General  Lunisolar 

preces.        precess.  100  m.         loo  ms.          100  n.  log.  100  n.        100 log.  loo  nt. 

"                "  "  s  "  s 

1725       502I.75  5036.22  4603.62  306.908  2006.l8  3.302369  133.745  2.I26278 

I75O       5022.30  5036.34  4604.32  306.955  2OO5.96  3.302323  I33.73I  2.126232 

'775     5°22-86  5°3647  4605.01  307.001  2005.75  3.302277  133.717  2.126186 

1800     50234i  5°36-59  4605.71  307048  2005.54  3.302231  133.703  2.126140 

1825     5023.97  5036.71  4606.41  307.094  2005.32  3.302185  133.689  -2.126094 

1850     502453  5036.84  4607.11  307.141  2005.11  3.302139  133.674  2.126048 

1875     5025.08  5036.96  4607.80  307.187  2004.90  3.302092  133.660  2. 1 26001 

1900     5025.64  5037.08  4608.50  307.234  2004.68  3.302046  133.646  2.125955 

1925     5026.19  5037.21  4609.20  307.280  2004.47  3.302000  133.632  2.125909 

1950     5026.75  5037  33  4609.90  307.327  2004.26  3.301954  133.617  2.125863 

1975     5027.31  503745  4610.60  3°7-373  2004.04  3.301908  133.603  2.1 2581 7 

2000     5027.86  5037.58  4611.29  307.420  2003.83  3.301862  133.589  2.125771 


Quantities  expressing  the  relative  positions  of  the  equator  and  equinox  at  any  two  epochs,  for 
use  in  the  trigonotnetric  reduction  of  mean  places  of  the  fixed  stars  from  one  epoch  to 
another. 

We  consider  the  quadrangle  formed  by  the  poles  of  the  equator  and  ecliptic  at 
two  different  epochs.  One  of  these  epochs  we  call  the  zero  epoch.  It  is  supposed  to 
be  that  for  which  the  position  of  a  star  is  given  and  to  which  the  various  quantities 
are  in  the  first  place  referred.  The  other  epoch  is  considered  to  be  a  variable  one. 
We  represent  the  respective  poles  of  the  ecliptic  and  equator  by  the  symbols  E  and 
P,  those  which  refer  to  the  zero  epoch  being  called  E„  and  P0. 

We  join  the  points  E0  and  P  so  as  to  divide  the  quadrangle  into  two  spherical 
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triangles  E„  P„  P  and  E0  P  E.  We  represent  the  parts  of  these  triangles  which  enter 
into  the  theory  as  follows: 

e0  =  E0  P0,  or  the  obliquity  of  the  zero  equator  to  the  zero  ecliptic. 

«!  =  E0  P,  or  the  obliquity  of  the  actual  equator  to  the  zero  ecliptic. 

{  rEP,  or  the  obliquity  of  the  actual  equator  to  the  actual  ecliptic. 

9  =  P  P,,,  the  distance  between  the  poles,  or  the  obliquity  of  the  actual  to  the 

zero  equator. 

^  =  Angle  P„  E0  P,  or  the  lunisolar  precession  upon  the  zero  ecliptic. 
9oo-£0  =  AngleE0P0P. 
oo°-£  =  Angle  E0  P  P„ 

A  =  Angle  E0  P  E,  the  planetary  precession  on  the  equator,  or  the  arc  of  the 
actual  equator  intercepted  between  the  two  ecliptics. 

*  =  ?  — A. 

The  formulae  embodying  the  application  of  the  preceding  quantities  to  the  reduc- 
tion of  the  fixed  stars  differ  from  those  of  Hessel,  used  by  Chauvknet,  in  that  the 
epoch  for  which  the  poles  E0  and  P0  is  supposed  to  be  given  is  the  same  for  which  the 
position  of  the  star  is  given.  In  the  other  theory  the  zero  epoch  is  a  constant  for  the 
whole  theory,  so  that  three  different  positions  of  each  pole  enter  into  the  theory. 
This  change  involves  a  slight  change  of  notation,  which  seems  advisable  in  order  to 
avoid  confusion.  The  notation  here  adopted  differs  from  that  in  the  concluding  pages 
of  "Elements  and  Constants"  only  in  that  I  have  here  put  £0  for  what  is  there  called 

The  reduction  in  question  involves  three  constants,  9,  and  z,  of  which  the  use 
is  as  follows.    We  put, 

«„,  £0>  the  coordinates  of  a  star  as  given  for  an  epoch  tv  which  I  call  the  zero  epoch. 
a,  8,  the  coordinates  as  reduced  to  an  epoch  t  by  precession  only. 

To  find  a  and  8  we  compute, 

a„  =  at  +  ^0 

p  =.  sin  9  (tan  8V  +  tan  -  9  cos  «o) 

.     ,        s        p  sin  a0 

tan  (a  —  Op)  =  —  — - — 

i  —p  cos  a0 

cos  -  (a  +  Oo) 
tan-(8  —  <$„)=      '  tanlfl 

2  I  2 

cos  -  (a  —  Op) 

a  —  a-\-z 

As  defined  by  these  equations,  the  constants  9,  £„,  and  z  are  functions  of  to  and  t, 
and  can  be  most  conveniently  expressed  as  functions  of  f0  and  t  —  f0)  because  then  they 
vary  very  slowly  with  t&    Their  values  for  any  pan  of  epochs  between  1600  and 
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2100  may  be  found  from  the  following  expressions,  which  are  given  for  three  values 
of  the  zero  epoch,     and  where  we  put 


250 

Zero  epoch. 

tt 

l600; 

e  — 

33  29  28.69 +  0.55 137*  — 0.12067* 

—  0.0002  T* 

1850; 

23  27  3 1.68 +  0.4078  —0.1207 

—  0.0002 

2100; 

2  3  2  5  34-56+0.2645  —0.1206 

//             tt  tt 

—  0.0002 

160O; 

I  2589.OI  T  —  6.66T*  —  0.027* 

1850; 

I2592.06    — 6.69     — 0.02 

2IOO; 

12595-»5  —6.71  —0.02 

//         //  // 

1600;  £+?„  = 

1 1 545-65»"  —  6. 1 27*  +  o.54rs 

1850; 

"SSMO  —6.14  +0.54 

2100; 

11556.98  -6.16  +0.54 

//          //  // 

160O;  f- 

45-3or  —  9-9lT*  —  °03  *3 

1850; 

33-54  -9-9'  -0-03 

2IOO; 

21.77  —9.92  —0.03 
ft           tf  tt 

1 60O; 

A  — 

45-30r  —  1483^  —  0.057* 

1850; 

33-54  —H-86  -0.05 

2IOO; 

21.77  —.14-88  —0.05 
//         //  // 

1600; 

&  = 

5750.177+  1.907*  +  o.28t» 

1850; 

5758.88  +1.89  +0.28 

2IOO; 

5767.61  + 1.88  +0.28 

tt         //  // 

1600; 

z  — 

5 750. 1  77  +  6.82 72  +  O.3OT8 

1850; 

5758.88    +6.83  +O.3O 

2IOO; 

5767.61     +6.84  +O.3O 

//               ft  tf 

1600; 

9  — 

5018.107—  2.667s  —  0.647s 

1850; 

5012.77  —2.66  —0.65 

2100; 

5007.44  —2.67  —0.65 

As  an  example  of  the  use  of  these  expressions,  let  us  suppose  the  epoch  for  which 
a  and  6  are  given  to  be  1875.  By  interpolation  of  coefficients  to  <o  =  1875,  we  find 
the  following  expression  for  £0,  z,  and  0 

tt  tt  tt 

1875;  Co  =  575975''+  1-89^  +  0.287s 
z  —  5759-751'  +  6.837*  +  0.30T8 
0  —  5012.247  —  2.66r*  —  0.657s 
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If  we  take  t=  1900  as  the  epoch  to  which  a„  and  S0  are  to  be  reduced,  we  should  put 

T  =  O.IO, 

which  would  give, 

// 

£0=  575-993 
s=  576.043 
9  —  501.196 

as  the  constants  for  reduction  from  1875  to  1900. 

A  check  on  these  values  is  afforded  by  computing  the  values  for  reduction  from 
1900  to  1875.    For  <o  =  1900  we  find  the  expressions 

//  //  // 

£0=5760.62  r-f  1.89  r2  +  o.28  T3 
z  =  5760.62  r  +  6.83  r2  +  0.30  T3 

9  ~  5O I  1 . 70  T  —  2.66  7*  —  O.65  T3 

For  t  =  1875  we  then  have  r  =  —  0.10,  and 

// 

£»  =  -  576.043 

*  =  - 575-994 
9  =  —  501.196 

If  we  put  the  values  of  £„,  z,  and  9  last  found  in  parentheses  we  should  have, 

(?.)  +  *  =  o 

0)  +  £0  =  O 
(9)+ 9  =  o 

It  will  always  be  advisable  to  test  the  values  of  these  quantities  in  this  way. 
The  same  test  may  be  applied  to  the  expressions  tabulated  on  p.  75,  by  taking  the 
epoch  t  as  1600,  1850,  and  2100  in  succession,  and  combining  each  with  the  other  two 
zero  epochs.  We  then  have  t  =  ±  1  or  ±  2,  as  the  case  may  be.  This  process  con- 
trols pretty  completely  the  entire  numerical  work  by  which  the  values  of  £<,,  z,  and 
9  on  p.  75  were  derived. 


